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I construct combined electric and magnetic field variables which independently represent energy 
flows in the forward and backward directions respectively, and use these to re-formulate Maxwell's 
equations. The emphasis is on detailed calculations, with a more general overview being published 
in Phys. Rev. A72 and arXiv. These directional variables enable us to not only judge the efi'ect and 
significance of backward-travelling field components, but also to discard them when appropriate. 
They thereby have the potential to simplify numerical simulations, leading to potential speed gains 
of up to 100% over standard FDTD or PSSD simulations. These field variables are also used to 
derive both envelope equations useful for narrow-band pulse propagation, and a second order wave 
equation. Alternative definitions are also presented, along with their associated wave equations. 

This report should be read along with the paper Phys. Rev. A72, 063807 (2005) "Theory of directional pulse propa- 
gation", by P. Kinsler, S.B.P. Radnor, G.H.C. New, for proper context. 

This document is primarily intended as a complete (as possible) record of the calculational steps that were necessarily 
abbreviated (or omitted) from that published work; it also contains a great many other relevant calculations along 
with some speculation. It is an edited version of a longer document from which on-going work has been excised; and, 
as a "work in progress", despite my best efforts, may contain occasional mistakes. Please contact me if you have any 
comments, corrections or queries. 
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I. INTRODUCTION 



We define field variables with directional properties, and use them to show how to optimise both our understanding 
of pulse propagation and numerical simulations. Simple plane-polarized versions of these were originally introduced 
by Fleck [Q, using e^^'^Ex ± n^^'^Hy. However, these were little used and the bulk of the paper primarily addressed 
generating ultra-short pulses by Q-switched lasers, discussing numerical techniques and simulation results. Here 
we derive wave equations and present results using generalised versions of these field variables. In the course of our 
investigations, we use these variables to illuminate several interesting features of z-propagated optical pulse simulations 
which are usually overlooked in a "forward-only" approximation. 

The main novelty of these new field variables is that they correspond to forward and backward directed energy 
fluxes; and are constructed from combinations of the electric and magnetic fields, along with the properties of the 
propagation medium. 

I also derive wave equations describing the propagation of these field components, in fully vectorised and plane- 
polarized versions. In addition to these exact re-expressions of Maxwell's equations, I also present second order and 
envelope- based wave equations for G^, in both stationary and moving frames. A further advantage of this approach 
is that it is as easy to include magneto-optic effects as electro-optics ones (i.e. dispersion and nonlinearity) . For 
example, I will later show that the plane polarized first order wave equations describing the z-propagation of in 
the spectral (cj) domain are - 

a.G± = T^toarf3r G± T * [G+ + G"] - '-^(3, * [G+ - G"] , (1.0.1) 

where ar = y^Ci, (3r = v^Mi contain the linear response of the medium, and ac contains nonlinearity that affects the 
electric field, and /3c is its magnetic counterpart. Using these forward and backward directed G^ variables has three 
main advantages: 

• The forward directed field variable G+ is the appropriate physical choice to use when studying the forward 
propagating part of the field - as is the case for most pulse propagation investigations. 

• Because we can choose to simulate only the forward directed field variable G+, we have half the calculations 
to do, and hence gain a speed advantage. By using a z-propagated PSSD algorithm^, we also get a fast and 
flexible treatment of dispersive and nonlinear effects; hence our simulations signiflcantly outperform standard FDTD 
approaches. 

• Most other approaches (even the recent |^) assume the backward propagating parts of the held are negligible. Since 
our G^ describes the backward directed part, we have a clear and physically appropriate basis on which to calculate 
it, and hence can (when appropriate) rigorously justify the common forward-only approximation. 

We can do even better by changing to a co-moving frame. This simplifles the G^ propagation, by removing the 
frequency-like oscillations, but has the price of making the G~ oscillations faster. However, in the limit where we 
decide to approximate by ignoring the backward propagating G^ part of the fleld, this price becomes irrelevant. 
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A. Why am I doing this? An anecdotal history 

Fleck's definitions were brought to my attention by Geoff New shortly after I started working with at Imperial 
in 2001. Because they gave rise to additional d/dt terms in the polarization, when compared to standard forms, he 
wondered whether there was a link with the Brabec and Krausz derivation[Q of similar-looking corrections to an 
envelope propagation equation. Although it turns out that this was not the case, I continued working on the idea. 
At first I focussed on generating a second order wave equation of , in analogy to the one usually presented for the 
E field. Although I generalized my GFEA derivation ^ to include Fleck- like terms, there seemed no practical use 
for the technique. 

In 2002 a PhD student in our group (JCA Tyrrell) started working on Maxwell's equations solvers, and applying 
them to pulse propagation in nonlinear media. He developed the PSSD method Q for solving Maxwell's equations, 
where the field is held as a function of time and propagated forward in space. This is common in envelope-based 
nonlinear optics simulations, but had not been applied to full Maxwell's solvers before. This got me thinking about 
representing Maxwell's equations using G^, and put the project on track to being of practical use, although there were 
many outstanding questions on interpretation of the fields and the wave equations While trying to clarify these issues, 
I became aware of work by Kolesik and co-workers ^), who projected out Maxwell's equations into a directional 
form. Their treatment motivated me on to vectorize the definitions, and get a PhD student (SBP Radnor) to start 
on simulations. As we worked on the paper eventually published as [^,1 further generalised the definitions to include 
the longitudinal parts of the fields, fully developed the variant forms of the variables, and derived envelope versions 
of the wave equations. 

As it stands, these kinds of directional wave equations, either using my G^ variables, or the Kolesik et.al. approach 
make traditional envelope theories based on the second order wave equation for E utterly redundant. This is true 
even for the generalised ones such as |^ which incorporate wideband corrections. Directional wave equations simple 
provide propagation equations without introducing the many approximations inherent in theories based on second 
order forms, as are detailed exaustively in j^, ^. This is treated in more detail in my report "Field and envelope 
methods in nonlinear optics". 

I am occasionally asked why I picked G^ instead of (e.g.) to name these directional fields - after all, F follows 
alphabetically from the electric field (letter) E^ and it is the first letter of Fleck's surname, so it might seem the more 
obvious choice. Unfortunately I can't remember why; but there was a reason why I passed over F , I think something 
to do with some other F I was using at the time. I should have picked if^, anyway. 



B. An alternative approach 

1. Kolesik, Moloney, and Mlejnek; PRL 89, 283902 (2002); and projection operators 



Kolesik et.al. |q] also use not-dissimilar combinations of field variables. They introduce projection operators for 



the forward and backward parts of the field in their cqn (KMM 1), which I rewrite slightly using u — kj 



D(k) =F sgn(fcz) u x H{k) 
bjyk) 



y^D{k) 

They then use these, with Maxwell's equations to generate their UPPE in cqn (KMM 7): 



Hik)) ± sgn(fcz)^^u X D{k) 
k 



dtDfik) = -iuj{k)Df{k) + -Lj{k)\PNL{D,k)^u[u- PNL{D,k) 

There are three main difference compared with my generalised Fleck-style G^ approach. They: 

1. describe things in terms of projection operators, 

2. project out both a directional D and a directional iJ, 

3. construct their (effective) directional variables in a different (but related) way. 



(1.2.1) 
(1.2.2) 

(1.2.3) 
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Another difference is that, when implimenting a solution, they keep the time derivatives on the LHS, and manipulate 
the V X terms on the RHS; so their UPPE gives a time derivative of the field, and so their fields are propagated through 
time. This contrasts with my approach which keeps the Vx on the LHS, and manipulates the time derivatives - 
hence my fields are propagated through space. Their choice is (most probably) why they choose to work with the 
displacement field D rather than E. In the time domain, the use of D avoids the need for the D — e* E convolution 
of the field with the permittivity, at least on the time-derivative side of the Maxwell's equations. They (presumably) 
remove the space derivatives resulting from the Vx terms when going into the spatial frequency (fc) domain, leading 
to the appearance of K terms on the RHS. 

However it is not so easy to avoid complication when solving Maxwell's equations. Their use of D simplifies the 
time derivative side, but there is a cost to pay on the space derivative (Vx) side, where they need to reconstruct E 
from D (see eqn (DMN 2)), and similarly evaluate the nonlinearity (see eqn (DMN 3)). Even then, just prior to eqn 
(DMN 7), Kolesik et.al. explain how they go into the frequency domain anyway! As an aside, they might easily have 
used B and not H in their equations, and so included the possibility of magnetic media in their model. 

In contrast, I leave the space derivative (Vx) sides alone, and work with the time derivative side. This has the 
advantage that when I go into the frequency domain, the convolution becomes a simple product, and I do not have to 
r esort to reconstructions of E and Pnl- It is interesting to compare the UPPE with my similar- looking eq n 



( 4.1.g ) t hat fo llows later, where they neglect to introduce the advantages of a co-moving frame (cf my eqn( 4.3.14 )) 



See also [II G, where I attempt a conversion between my formulation and theirs. 

SUMMARY: Maxwell-like theory with combined G-like variables, but no co-moving frame. Temporal propagation 
only. 



2. Kolesik and Moloney; PRE 70, 036604 (2004). 

This work 1^ is an extension of the method proposed in and contains a variety of applications. It discusses 
both t and z propagated approaches, and includes decomposition into transverse modes. They develop an envelope 
version of their wave equation, but do not discuss using multiple envelopes. In general they only solves for E, not 
both E and H separately, there is also no examination (or calculation) of the significance of backwar d pro pagating 



fields. They separate the polarization P from displacement D; and includes current terms. See also [II G, where I 
attempt a conversion between my formulation and theirs. 



C. Moving frames 

In Fidel et.al. 1997 ||lo|| "Hybrid Ray FDTD moving window approach to pulse propagation"^ they look at a variety 
of moving window and moving frame approaches to solving for pulse propagation with Maxwell's equations using 
FDTD. They do introduce a co-moving frame, and derive the appropriate discretized Maxwell's equations; but they 
do not use combined E and B variables as done by Fleck, and as done in the following sections here. 
SUMMARY: Maxwell-like theory with co-moving frame, but no combined G-like variables. 



D. Beltrami Fields, Photon Wave Functions: E + iB 



The fields I constuct here bear a superficial resemblance to Beltrami variables 1 11 , ij, |T3, 14 , which are defined 
along the lines of Q = ^/eE + i^JJlH . However, they differ in two important respects. First, for a given Beltrami Q you 
immediately know what the component E and H fields are, whereas you require both G^ to do the same. Secondly, 
Q does not assume any preferred direction, whereas the construction of G^ requires a direction u to be chosen. It's 
interesting to wonder whether the E ± iH approaches might generalise to something like E ± e^'^H). However, since 
E ± iH is unique for a given E and H; and E zL H is not unique without both components, so it's not obvious how 
it might work. 

Note that some discussions of Beltrami fields refer to "force free" systems, where V x _B = aB, which are not (as 
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far as I can see) applicable to an optics regime. 
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II. THE MATERIAL PROPERTIES: REFERENCE AND CORRECTION PARAMETERS 



Maxwell's equations describe propagation of the EM field through a medium described by permittivity e and 
permeability /i. In addition, it is common to add in a "polarization" term P, which is used to relate the displacement 
field D to E in the usual way: D = e * E + P. In my treatment, I hide this polarization inside the definition of e 
(writing instead just D = E for the same system); similarly any magnetic polarization would be hidden inside /i 
and not written out separately. Note that I do make a split of the total permittivity (and also the permeability) into 
two pieces, a "reference" part and a "correction" part, and these will not necessarily be the same as the traditional 
division. 

Note 1: I use the "*" notation for convolutions a{t) * 6(t) = J a{t — t')h{t')dt' . 
Note 2: Quantities in the frequency domain are indicated by tildes. 



A. Reference and correction contributions 



We first denote the linear response of the medium to be e^, and its nonlinear response to be e^L- The nonlinear 
term may include some time response, such as e.g. €NL{t) = / e3{T)E{t — T^dr. The effect of the total permittivity 
can be expressed in either the time domain or the frequency domain (indicated by tildes). It is 

e{t)^E{t) = tL{t)*E{t)+eNL{t)E{t) (2.1.1) 
or i[uj)*E[uj) = iL{i^)E[uj) + lNL{i^)*E{uj). (2.1.2) 

Similar expressions can be written down for the case of magnetic dispersion or magnetic nonlinearities, 

fi{t)*E{t) = fiL{t) * E{t) + ^lNL{t)E{t) (2.1.3) 
or fl{uj)*E{uj) = jjLL{'^)E{uj) + jlNhi'^) * E{uj). (2-1.4) 

The definitions of (and their generalized vector counterparts , introduced below) depend on the properties 
of the propagation medium through the permittivity e and permeability /i. In principle it might seem attractive to 
define G"^ using the exact values oi e,ii (including the nonlinearity) , but we will usually want to be able to uniquely 
reconstruct the fields E and H from our new fields. We do this by including as much as possible of the dispersive 
properties in the reference parameters (and hence the definitions of G^), leaving nonlinear properties (and potentially 
some residual dispersion) for correction terms. For example, if the true dispersion was multi-valued or zero in places, 
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it is possible to define the reference to include an approximation to the true dispersion without those inconvenient 
properties, and use the correction term to compensate. 

The frequency domain is the starting point for calculating the necessary parameters. This is because the formulae 
are simplest in the ixi domain, where any linear response is described using products rather than convolutions. This 
will mean that we calculate the time-domain parameters by (back) Fourier transforming the frequency domain ones. 
The frequency domain definitions are 

^ er{uj)E{uj) + lc{uj)*E (2.1.5) 
= al{Lu)E{uj) + &r{Lo) aduj) ^ E{uj), (2.1.6) 
jl*H = jlr{uj)H{uj)+jlc{io)*H{uj) (2.1.7) 
= pl{uj)H{uo)+pr{uo) h{^)*H{u;), (2.1.8) 

where the correction parameters Cc and flc represent the discrepancy between the true values and the reference. 
The smaller these correction terms are, the better the match. By using these frequency dependent parameters in 
the generalized definitions of G^, we will be able to propagate pulses using only the G'^ variable, a gain in both 
mathematical simplicity and computational speed. Note that since the definitions of depend (roughly speaking) 
on the square roots of e and /2, we have introduced the a and [3 parameters, which will feature prominently (along 
with their time domain counterparts a, /3), in the generalized definitions of G* that follow. 
To recap, I split the medium properties into two parts: 

First, there are the "reference" parts and /i^. These are the contributions from the material properties which 
will be incorporated into the definition(s) of G^. These will only contain linear response terms, to guarantee we can 
reconstruct E, H from G* 

Second, there are the "correction" parts ec(t) and fj,c{t), which will contain the difference between the chosen 
reference and the true behaviour of the material. These will only contain any nonlinear response terms, but may also 
contain an residual linear response also. 

We might also choose to split the correction terms up into more pieces - e.g. a dispersive part, and a nonlinear 
part: e.g., 

ac*E = a^E + a^^ * E. (2.1.9) 
B. Nonlinear ity 

Since it is usually impractical to include nonlinearities in the reference parameters, these will normally appear 
in the correction terms ec, Hc- As an example, consider a n-th order (electric) nonlinearity, in which case f-dt) = 
x("Hi) * -B(t)""\ and 

ar{u;)a,{uj) = :r [x^^HO * ^W""'] (2.2.1) 

cie(^) = [ar{iu)]-\j [x^"Ht)*E{tr-'], (2.2.2) 
ac{t) = {[ar{Lo)]-\x^'^\uj).j[E{tr-^]} , (2.2.3) 
where (see later) E{t) = i {S"-^ [a~^] } * [G+ + G"] (2.2.4) 



where ?[...] is the Fourier transform (FT) from time to frequency, and E{t) can be found from eqn. (3.3.3). If 
the reference parameters a,, contain dispersion (which will be the typical case), we can see from eqn. ( 2.2.2| ) 
that this will make adw) dispersive even if x^"^ is instantaneous. In the case of an instantaneous nonlinear- 
ity, this adds more computational work (an extra two FTs), although for non-instantaneous ones we needed the 
FTs anyway. If the nonlinearity is instantaneous and the reference parameters are non-dispersive, we have simply 
^i(t) = a-".x(").2-"+i [G+ + G-]""\ 



In practise, therefore, calculation of an instantaneous nonlinear term will involve calculating x(")i?(<)", fourier 
tranforming (denoted by J) the result, then dividing it by the (possibly frequency dependent) reference parameter 
Ur- If the reference parameter is a constant, then the fourier transform is redundant and the nonlinear step can 
be handled while remaining entirely in the time domain. Note that there is no way we can avoid convolutions if 
the reference ar contains any frequency dependence - even if our nonlinearity is instantaneous. This is because to 
calculate ac I have to remove the reference ar from ic] since adw) = ec{uj)/ar{uj). 
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III. DEFINITIONS OF G± 



Here I go significantly further than Fleck[0 in two respects: I vectorise the definitions, and allow for dispersive 
effects. An important addition is the inclusion of a longitudinal field component, which is required in order to retain 
a full vector description of the EM field. Originally I tried to generate a motivation for the directional fields 
described above by taking Maxwells equations and moving to a co-moving frame. However, it seems now that while 
they certainly seem to describe parts of the field propagating in opposite directions, they do not seem tied to any 
particular choice of frame. 

After a brief discussion of what Fleck did, I define a vector form of the Fleck-style fields G*. This requires the 
definition of a "direction of propagation" denoted by the constant unit vector u. I can then achieve the swapping of 
the transverse components of H (since ~ ± Hy) by using the cross product of u with the magnetic field variable 
H. 

Secondly I use "reference" material parameters which I use to construct my G* fields. This allows me to distinguish 
between the parameters used in my construction of G^ from the actual material properties, which may be too 
complicated to use in that context. Generally, I include all of the linear response of the medium in these reference 
parameters, which makes it easiest to define G^ in the frequency domain. This is because the linear time-response 
convolution in the time-domain becomes a simple product in the frequency domain. A good (accurate) choice of 
reference medium for a particular problem will ensure that the contribution due to the backward directed part of 
the field (G~) is small. In general we would want the reference part to include as much about the material as is 
practicable, to ensure the smallest possible correction term(s) for propagation. 

Note that whether or not the e, fj, contain dispersion, we here restrict ourself to the case where they are scalar. Also, 
I let the argument of the vector field quantities tell us whether we are in the time or frequency domain, to avoid too 
much notational clutter. As a final note, in the following calculations the assumption will be made that the media 
are bulk, without (implicit) interfaces hidden in e or /i. This means that dze = dzH = is assumed. 



A. Fleck's approach 



Fleck defines a relative permittivity e, and magnetic permeability in a dispersionless host medium - it is clear 

1/2 

that these are relative parameters because he later defines 77 = (e/z) as the refractive index of the of host medium 
(just after eqn (F1.4b)). He allows for a Conductivity cr; and a Polarization P which characterises the active atoms 
and which may be amplifying or absorbing. The electric field is Ex and the magnetic field By — iiHy. As usual, the 
speed of light c is just related to the permittivity and permeability of the vacuum, i.e. (cqMo) ■ Fleck simplifies 
plane polarised Maxwell field equations, treating them only along the z direction, resulting in his eqn.(Fl.la,b) 

dt-Ex = -dzHy~—dtP^—aEx, (3.1.1) 
c c c 

dt!^Hy = -dzEx. (3.1.2) 

Then Fleck defines new combined EM field variables E^ in eqn.(F1.3), but I write them G^ to avoid confusion 
with the E field variables. 



E^ 



G^ 



(3.1.3) 



Because of the relative phases of E and H in these definitions, if the the two fields G+ and G^ are constants 
multiplied by a standard exponentially oscillating carrier wave, they are associated with an energy flux (i.e. Poynting 
vector) in opposite directions along the z axis; this will presumabl y also hold nearly true when replacing the constant 



amplitude with "slowly varying" pulse envelopes. See also section VIII 



Fleck uses his definitions to obtain equations of motion for his field variables G+ and G {E ), recorded in his 



eqn.(F1.4). They are 



^dtG+ 
c 

'-d,G- 



dzG+ 
dzG- 



Anil 



1/2 



-dtP 



1/2 



-dtP- 



c 
2tt(j 
c 



(?) 
(f) 



1/2 



1/2 



(G+ + G" 
(G+ + G- 



(3.1.4) 
(3.1.5) 
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I now scale eqns.(|T|,|T|) (Fleck (F1.4a, F1.4b)), setting G"± = {ri/c)G^, P' = {Atth^^^ /c) P, and a' 



1/2 

{2'Ka/e) . In the following, I will drop the prime (') marks and use G for G' (etc) in order to reduce the visual 
clutter and make equations easier to read. For this subsection only, from this point onwards each occurence of G^, 
P, or a should be understood to refer to the prime variables G'^, P' & cr' - unless explicitly stated otherwise. The 
scaled equations are 



5tG^ 
dtG- 



9.G^ 
dzG- 



-dtP 
-dtP 



a{G+ 
a{G+ 



G- 
G- 



(3.1.6) 
(3.1.7) 



B. Alternative forms 

The "standard" definition for directional fields handles the electric field E most easily; although rather non- 
intuitively we will find that it makes it hard to include longitudinal E components. Fortunately, other definitions 



are possible , an d thes e follow the standard "Primary definition presented next in subsection III C , in subsections 

iiiDl fuM and fnH. 



Use of an alternative definition can protentially shift the non-transverse effects to a field variable where it causes 
less trouble. Note that circularly polarized forms can be generated for all types the same way as is described for 



primary- i? type G^ variables in subsubsection III C 2 



These alternative expressions are represented only in the frequency domain; time domain forms are presented for 
primary-iJ only; for the other forms it is a simple matter to construct their time domain representation. In the case 
of dispersionless reference media (i.e. constant 6^(0;), /i(tj)), the time and frequency forms look very similar. The set 
of alternative definitions was (partly) inspired by PDD suggesting (private communication, 24 February 2003) trying 
alternative combinations - such as swapping D for E and/or B for H. 



C. Primary E: the standard vector form 



This vector definition of directional fields most closely mirrors the simple form suggested by Fleck. For completeness, 
I will write them down in both time and frequency domains. Note that the use oi u x H means that the G^ will not 
contain any information about u ■ H, the longitudinal part oi H ~ thus we insist G^ be magnetically transverse and 
define a G° . The vector fields G^ are (time domain) 



so 



G^{t) 
ar{t) * E{t) 



ar{t) * E{t) + u X (ir{t) * H{t), 

U- I3r{t)* H{t) 



and 



UXPr{t)*H{t) = - 



G+(i)+G-(t) 
G+(0-G-(t) 



or, in the frequency domain, where we can avoid the convolutions 



and 



so 



or 



G±(c.) 

ar{uj)E{Lu) 
u X [3r{uo)H{uj) 



ar{l^) E{uj) ± U X \/ /i,-(ti') H{llj), 
U ■ iSriuj) H{uj) , 



1 

2 

1 

2 

1 ^ 
-u X 
2 



G+{uj) + G-{uj) 
G+ [lu) - G- {uo) 
G+{uj)-G-{u) 



uG°, 



(3.3.1) 
(3.3.2) 

(3.3.3) 

(3.3.4) 

(3.3.5) 
(3.3.6) 

(3.3.7) 

(3.3.8) 

(3.3.9) 



since uxux A — [u-A] — A. In general usage I will remain in the frequency domain, because there I avoid complications 
due to convolutions (or worse, deconvolutions); in addition the notation is simpler. 

Definition: THF - transverse H (magnetic intensity) field (i.e. u - H = 0). 



A derivation of the wave equations using this form is in section IV 
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1. Divergence 



In the normal description of EM, the divergence V ■ D — p. However the picture here is based on E, so we need 
a way to construct D from E. We follow the following iterative procedure, which should work for the usual case of 
weak nonlinearities 



D{uj) = alE{Lu) + OrCtc * E{lu) 
alE{uj) 



D{lu) — OLrOLc * E{lj) 



E 



(3.3.10) 
(3.3.11) 
(3.3.12) 
(3.3.13) 
(3.3.14) 
(3.3.15) 



c>L^ Uc* D{lj) + OLc* otc* D{uj) — etc* etc* OLc* E{uj) (3.3.16) 



(3.3.17) 



Thus, 



OO 



v-i?(w) = y ■^a-;'^-''[-ac*X D{uj) 

OO 

= ^5;2-»[_5^*]*V-^(cc>) 

i=0 

OO 

= E [-"c*]* P{UJ) 



with the shorthand notation 



i=0 



(3.3.18) 

(3.3.19) 

(3.3.20) 
(3.3.21) 



It is worth noting that in the alternative constructions following this "Primary i?" one, it is far simpler to calculate 
the divergence. In any case, we can now write the divergence of as follows 



~-2 ~ 



V • Z3 ± PrH • (V X u) T PrU -(v xH^ 



'lUja^E + J 



a^arP -ii —a'^'a^^ PrU ■ 



=F PrU ■ J. 



(3.3.22) 
(3.3.23) 

(3.3.24) 
(3.3.25) 
(3.3.26) 



Thus I never use the (no monopoles condition) V • i? = 0. Clearly, if are pure transverse (not just THE, which 
is built into the definition, but no longitudinal E fields either), then 



NB: 



V • G^(ijj) = a^^&rP 



V • G+(cj) - V • G"(a;) = a^^^rp + ya^a^^/S^w • [G+(cj) + G"(w)J - • J(w) 



~-2 - ~2~-l a ^ 

a-^ ttrP + Pr'fJ- 



G+ + G- 



— PrU ■ J 



+iuja-^dr PrU- G (lu) + G (lu) 2PrU ■ J{lij) 



(3.3.27) 

(3.3.28) 
(3.3.29) 
(3.3.30) 
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2. Circularly polarized 

Handling circularly polarized light is easy using fields. Just as for E and H fields, you define 

E = xEx cos{u!t —(/))+ yEy cos{ujt), 
H = —xHx cos{ujt) + yHy cos{ujt — cf)) 

we get 

E±zxH ^ x[Ex±Hy] cos{ujt) + y[Ey ± H.^] cos{ujt - 0) 
= xG^ cos{ujt yGy cos{ujt) 

Similar definitions will hold for choices other than the standard "primary £"' G^ fields. 



(3.3.31) 
(3.3.32) 

(3.3.33) 
(3.3.34) 



3. Comparison with E and H 



Assume 



so 



G+ 
E = 

E/H = 



A 



G- =Q 



t{kz—(jjt) 



2ar 

A 2(3r 

2ar' A 



2/3. 



l{kz—LL!t) 



Compare with the results we get from an E carrier-based approach 

E = Be'('='=-"*\ 
to get H we use /ii9t-ff ~ d^E 

so ^l^tH = ik.Be<'''-'^*\ 



so 



E/H = bJ ^ 



/3 



olB 



a 



(3.3.35) 
(3.3.36) 

(3.3.37) 



(3.3.38) 
(3.3.39) 
(3.3.40) 

(3.3.41) 
(3.3.42) 
(3.3.43) 



D. Primary H 



If we wanted to swap the roles of H and E, we might instead define the field variables as 



so 



or 



and 



G'±(c.) 

G'°iLj) 

u X ar{uj)E{ijj) 

ar{uj)E{LLl) 



u ■ (5r(w) E{lo), 
G'+{uj)-G'-{lu) 

G'+{lu)-G'-{uj)\ +uG'° 



1 

2 

1 ^ 

-u X 
2 

1 
2 



G'+(cj) + G'-(cc;) 



(3.4.1) 
(3.4.2) 

(3.4.3) 
(3.4.4) 
(3.4.5) 



We would do this so that the magnetic terms not included in the reference medium (i.e. in the definitions of G'^) 
would lose the non-transverse correction terms (see later) , although terms like those would then appear on the electric 
terms (and I would define a transverse electric fields (TEF) approximation instead). 

A derivation of the wave equations using this form is in section 
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Divergence 



Here we need to calculate B from iJ, giving us something like B = f3-^'^H; just as for the Primary E form where 
we needed to calculate D from E. However, rather than repeat an almost identical calculation, I simply note that 
V • i? = (no magnetic monopoles), which will cause the term to vanish anyway. 

The divergence calculation is 



P^^Pr^ ■ B ± ■ {ux E 



T aril ■ (v X Ej ± ctrE (V x u) 



G+ [lo) + G- {lo) 



(3.4.6) 
(3.4.7) 
(3.4.8) 
(3.4.9) 
(3.4.10) 
(3.4.11) 



So no longitudinal magnetic intensity field H means zero divergence. Clearly, if G'^ are pure transverse (not just 
TEF, which is built into the definition, but no longitudinal H fields either), then the divergence is zero. 



E. Primary D 



If we were more interested in D than E (or B than H), we might instead define variables as 

F±(w) = a-^{uj) D{uj)±ux i3;^{uj) B{uj), 
F°{cj) = u- Biu), 



SO D{uj) 
and ux B{uj) 
and B[ijj) 



2 
2 



-U X 



F+{lu) + F-{lu) 
F+(w) - F-{u;) 

F+{uj) - F-{uj) + l3r{Lo)uF°{uj). 



(3.5.1) 
(3.5.2) 

(3.5.3) 
(3.5.4) 
(3.5.5) 



We would do this so that the displacement field (and magnetic fields) could be easily reconstructed from F^; and, 
in analogy to the THE and TEF mentioned above, we would use a transverse magnetic B field approximation (TBF). 



A derivation of the wave equations using this form is in section VI 



1. Divergence 



V-F^{lu) = V-a-^D{Lu) ±V -ux f3-^B{uj) 

= a-^V-D ±(3~^B-{\7 xu) T /3.7^w • f V x B 



a-'p ±0 T Pr P u - [\/ X H 



p ± —OtrPr P U- 



F+{iu)+F-{Lo) TPr'P^Ji^) 



(3.5.6) 
(3.5.7) 

(3.5.8) 

(3.5.9) 

(3.5.10) 
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Thus I never use the (no monopoles condition) V • i? = 0. Clearly, if are pure transverse (not just TBF, which 
is built into the definition, but no longitudinal D fields either), then the divergence is depends only on the charge 
density 



NB: 



V-F+-V-F- = iLoarPr^P'^u- F+ + F' -2(3-^0^. J 



(3.5.11) 



(3.5.12) 



F. Primary B 



If we wanted to swap the roles of B and Z?, we might instead define variables as 

F'^{uj) = p;^{uj) B{oj)±uyi D{uj), 
F'°{uj) = u- a-^{Lo) D{lu), 



SO 



ux D{uj) 



or 



and 



D{lo) — — ^r-^u X 



Biuj) 



2 



F'+ {uj) ~ F- (uj) 

F'+{uj) - F'{uj)] + ar{uj)uF'°{uj) 
F'+{lu)+F'-{lu) . 



(3.6.1) 
(3.6.2) 

(3.6.3) 
(3.6.4) 
(3.6.5) 



In analogy to the THF and TEF mentioned above, we would use a transverse displacement field approximation 
(TDF). 



A derivation of the wave equations using this form is (not yet) in section VII 



1 . Divergence 

V-F'^(cj) = V ■ i3-^B{uj)±V ■ (ux &-^D{uj)] 

= T a^a^^M ■ (\/ X ± ^^a-^E (V x u) 
— ^a^a^^u- {lojB^ 



F+{uj) + F-{uj) 



(3.6.6) 
(3.6.7) 
(3.6.8) 
(3.6.9) 
(3.6.10) 



Clearly, if F'^ are pure transverse (not just TDF, which is built into the definition, but no longitudinal B fields 
either), then the divergence is zero. 



G. Comparison to Kolesik et.al 



(needs updating with alteration of F, F' definitions at 20050727/8) 
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Now let us compare the equations Kolesik equations for D with those for the fields. They have 



D{k))^ 



D{k)) =F sgn(fcz) u x H{k) 



u{k) 



1 



c 



= ^/JID =F fI^/eu x H 
= y/JlD T V^u X B 



(3.7.1) 

(3.7.2) 

(3.7.3) 
(3.7.4) 
(3.7.5) 



Barring my dropping of the sgn(fcz) factor, their can clearly be identified with my "primary D" formulation 
of directional fields (see III E| ). Of course, there are differences, notably that their dispersion is only applied to the 
magnetic field (either H or B); my "primary i?" formulation applied e dispersion to E only, and /i dispersion to H 
only. 

They also present an equation for a directional H which corresponds to the "primary _B" form {F') as opposed to 
the D (F) one. 



H. Moving Frame 

It is often useful to transform the following wave equation derivations into a moving frame. The frame will be 
defined by the speed c/ = l/a//?/ in analogy to the reference and correction parameters above. A sensible choice of 
c/ is probably helpful, e.g. Cf = Cr or possibly c/ = Vg, the pulse group velocity. 

The (non relativistic) z-direction frame translation itself is 

t' = t-z/cf (3.8.1) 
z' = z (3.8.2) 



so that 

dv - dt (3.8.3) 

d,, = -cj'dt + d. (3.8.4) 

When doing vector calculations, I will need a vector form of the frame translation. I use 

V xQ V xQ-afPfUxdtQ (3.8.5) 

Assuming u is along the z-direction - 

VxQ ^ X {dyQ, - d.Qy) + y (9,Q, - 9,Q,) + z {d^Qy - dyQ,) (3.8.6) 

= X {dyQz - dz'Qy ~ afPfdtQy) + y {dz'Qx - dxQz + aff3fdtQx) + z {dxQy - dyQx) (3.8.7) 

= x {dyQz - dz'Qy) + y {dz'Qx ~ dxQz) + z {dxQy - dyQx) - xaf(3fdtQy + yaf(3fdtQx (3.8.8) 

= V xQ-xafl3fdtQy + yafl3fdtQx (3.8.9) 

I will often find it useful to define the ratio of the reference and frame speeds: 

^ = UfPf/UrPr (3.8.10) 

This moving frame is very useful in a space-propagated model since the pulse is held as a function of time, and so 
it will stay (nearly) centered while propagating forward. Although it has no sensible limit as the frame speed tends to 
zero, we can recover that case by setting a//?/ = (or ^ = 0) and rewriting the primed variables as unprimed ones. 
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1. Frequency domain 



In the frequency domain, the space derivative transforms like 

d^' = iLj&f(3f+d^ (3.8.11) 

which seems to leave open the intriguing possibility of a frequency dependent ( "dispersive" ) moving frame by 
allowing af and /?/ a frequency dependence. I could put all the linear evolution in the reference, and match it to the 
frame (so that af = ar and Pf = (3r and ^ = 1). This means that the only alterations to the pulse profile would be 
those due to nonlinear effects. 

We might hope to use a dispersive frame in the following way. Imagine a pulse being propagated in a strongly 
dispersive medium with some nonlinearity; and that the pulse might be expected to broaden by a large factor of 
(e.g.) twenty during propagation. This would mean we would need a much wider time window in our simulation than 
suggested by the initial pulse width. In a dispersive frame matched to a dispersive reference medium, our simulated 
pulse would not broaden at all due to dispersion, so we then hope retain a narrow time window throughout the 
simulation, hence saving computer time. The drawback is that our nonlinear term absorbs the reference dispersion, 
and this may cause significant pulse broadening of itself, perhaps of similar extent to the dispersion would. Assuming 
it did give an advantage, at the end of the simulation we would have to take our dispersive-frame pulse and convert 
back to a non-dispersive frame to see the "lab" pulse profile. Fortunately the computational cost of this conversion is 
minimal. 



I. Pseudo-Beltrami forms 

Weiglhofer and Lakhtakia present a number of quantities related to the EM field using their Beltrami style 
definitions. If I rearrange the factors of /i's and e, I might rewrite these as 

free field: e^'^Q± = e^'^E ±i^i^'^H (3.9.1) 

coupled field: e^/^F± = e^^/^D ± z/i^^/^^ (3.9.2) 

currents: ^l-^/'^W± = -\r^l'^ J^±it-^l'^ Je (3.9.3) 

charges: /i~^/^w± = -^i^^^Vm ± «e"^^Ve (3.9.4) 

Analogously, and since we use ux the curl-i/ equation (which has Z?, Je on the RHS) - 

free field: G± = e^^^E ± ^^^"^u x H (3.9.5) 

coupled field: F± = e^^/^u x i) ± ^"^/^^ (3.9.6) 

currents: K± = ^^/^m x X ± e^/^ J„ (3.9.7) 

charges: L± = e^^/^^e ± A^^^^Vm ?? (3.9.8) 

It might also be interesting to construct Beltrami-like combinations of the directional fields, i.e. 

g'=^ = G^+iG^. (3.9.9) 

J. Miscellaneous: Force on a moving charge 

If I define a velocity v — vu, there is a curious similarity between the force on a moving charge and the construction 
of my field variable. We have - 

-F = E + vx B (3.10.1) 
Q 

gl/2 

F ^ e^^^E + ve^/^ux fiH (3.10.2) 

q 

= + i,ei/ V^'^ X (3.10.3) 

= e^/'^E +^ux n^/^H (3.10.4) 
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So if t; = c, 



cl/2 



= G+ 



(3.10.5) 
(3.10.6) 



This seems to be a remarkable coincidence! 
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IV. FIRST ORDER EVOLUTION EQUATIONS: PRIMARY E 



Here I derive coupled wave eq uatio ns for the fields defined in III C , I also discuss how things look in various frames 



of reference, and in subsection IVB generate first-order differential equations for forward and backward propagating 
fields. In the following I do not introduce any envelope-carrier decomposition for the field variables; alt houg h the 
introduction of a moving frame plays a related role. Note that the Vector Stationary Frame subsection ( IV B ) also 
contains the plane-polarized special case, as well as a transformation to the moving frame as well; and should be 
considered to contain the authoritative results. Note there are still some sign discrepancies between the vector 
calculations and the plane polarized ones. 



A. Plane-polarized stationary frame (Example only) 



Note that is better to follow the (authoritative) vector stationary frame derviation (see IVB ) instead. However, 
his is a simpler derivation that gives a nice overview of the approach. It is superceeded by the full vectorized (and 
moving-frame) derivation(s) below. I begin with the source-free Maxwell's equations - 



dz 
dExjt) 
dz 



(4.1.1) 

(4.1.2) 



Because the convolutions e^Ej. and fJ.*Hy complicate the analysis, I will fourier transform into the frequency domain 
where their linear parts become simple products. Using the definitions in subsection |l| and carelessly dropping the 
tildes off the a and /? parameters, I can proceed with the calculation ~ 



oz 

-^PrHy 

az 
d 

-TT-arEx 

OZ 



— iLoeEx 




ibj [alEx + 


arCic * 


— ILOjlHy 




= lUJ [plHy + 




= lUJarPr 


UrEx 


-I- ttoPrCtc * 


arEx 


= lUOarPr 


PrHy 


+ VjjarPc * 


PrEIy 



(4.1.3) 
(4.1.4) 
(4.1.5) 
(4.1.6) 



Adding and subtracting eqns. ( 4.1.5 , 4.1.6 ), and using the fourier transform of the definitions for = arExi(3rE[y, 
we get 



— — -G^ = ±lUjar(3rG^ + ILOarPc * f3rHv ± lL0(3rCtc * CtrEx 

oz' 



^Y/3c*(g+ 



G- 



Pr ( 



G+ + G- 



(4.1.7) 
(4.1.8) 



The "non-magnetic" limit (strictly, no magnetic corrections to the frame), where /3c = is 

^G± = ±iujarl3rG^ ± iw^ac * f G+ + G" 
oz' 2 V 



(4.1.9) 



T he refer ence-only limit of the above equations (with ai ~ 0; /3i = 0) is exactly as you would expect (see 
eqn.(^l.l.lOD) - 



ztlUjarPrG 



(4.1.10) 



Clearly, even if ac <C a^, both G+ and G evolve on simila r scales, as specified by the leading RHS term. It is 
interesting to compare Kolesik et.al.'s eqn(KMM 7) (i.e. (1.2.3)) with my eqn (4.1.8) above. Of course their equation 
has the G^ terms projected out, but otherwise the two are similar in form. 
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B. Vector stationary frame (Authoritative) 

Here I avoid calculations using the x, y, and z components of the field, and retain a fully vector description. Using 
the definitions in subsection I can proceed with the calculation - 



into temporal-frequency space 
swap transverse components — 

premultiply — s- 



V X H{t) 

V X E{t) 

V X H{uj) 

V X E{ijj) 
ux (v X H 



sum-and-difference 



V X 

ux (V X Prltj 

V X a^E 
V X arE ± ux (v X l3rH 



+dte * E[t) + J, 

-dui*H{t) (4.2.1) 

{ujf * E{uj) + J(cj), 

+ILO Piuf ^ H{uj) (4.2.2) 
—luj * {u X E^ + u X J, 

+ILJ 0^ *H (4.2.3) 

— ILO (irCp' * X itj + UX PrJ, 

+ILJ arfi^ * H (4.2.4) 

+IUJ arf3'^ * H ^ lUJ PrCt^ * {u X Ej 

± uxf3rJ, (4.2.5) 



The vector G fields (defined in eqn (3.3.1 )) are 



arE + U X PrH 



(4.2.6) 



This means I need to convert both the second term on the LHS of the sum-and-difFerence equation above, as well 
as the RHS. It is most imortant for the LHS to be simple, because this will define the type of propagation specified 
by the RHS. To convert the LHS, we need two vector identities, as used in subsection XH to get eqn.(12.0.1C). So 
using 



u X V X - V 



[u-h) 



= V X [ux H 



(4.2.7) 



If I retain the u ■ PrH term , which will be zero for transverse H fields (THE), then 



+iuarl3 *H =F lujPrCt *(uxE] ± UX(3rJ 



V X arE ± u X (^V X i3rH 
V X arE ± V X f M X l3rH) ± V f u • PrH 



-\-lL0arl3 * H =P ILoPrCt * 



[ux E^ 



± ux PrJ 



(4.2.8) 
(4.2.9) 



V X 



arE ± {ux PrH 
V X G± 



-lujarfj'^ *H T lioPrO^ * [ux E^^V [u- prH^ ± w X /3r(a.2.10) 
j^^arf3^*H T Pra^ * [ux E^^ [u- f3rH^ ± m x /3r/4.2.11) 



Now note that eqn.( 12.0.12| ) means that 



u X 



and so 



ux H 
uxG^ 



u-H 



u-H 



u X arE ± u X 



u X (3rH 



u X arE =F (3rH ± 



U-(3rH 



(4.2.12) 
(4.2.13) 
(4.2.14) 



So, 

V X G± 



U-H 



ar0 * [ux 



\^[3ra^ *[uxE^± 



ux H 
ux H 



T Pra^ [ux Ey^^V [il- PrH^ ± u x (3rJ (4.2.15) 



)} 



+ lUjarP * 



U-H 



U ^v{u-l3rH] ± zi X /3r44.2.16) 
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I now separate the interaction parts (depending on ofc, Pc) from the vacuum-like parts (depending on g ^i 0r) i and 
t hen s ubstitute (as far as possible) expressions containing G rather than E or H, by referring to eqns ( 3.3.3 ) and 
( |3.3.4| ). Hence 



V X G"^ — ^lUJUrPr |u X arE ±U X 



U X 



=F VjoPrUc * X UrE^ — lUOaj-Pc * ( M X 



TlUJarPr U X G + lUCirfirU G° =F VG° 

ILUar 



lUjj3r 

T — ;:; — oic * I w x 



G+ + G- 



U X PrH 



f3c* (ux 



U ■ (3rH 



u ± uxPrJ (4.2.17) 



+ iujarl3c*{uG°) ± w X /3r^4.2.18) 



Here the transverse (THE) and longitudinal parts decouple, since is guaranteed magnetically transverse, and 
uG'° is magnetically longitudinal. The two decoupled equations are - 



V X G^ ~ ^luarPr u X G"^ q= 



lLjj3r 



-ar * \ u X 



G+ + G- 



G+ -G- 



±VG° = +lUJar(irUG° + lUJUrfic* {u G°) . 



± ux(3rJ (4.2.19) 
(4.2.20) 



We see from the following tV B 1 that these two equations are sufficient to describe the fields as long as there is no 
lo ngitudin al E component (i.e. G^ is transverse) or longitudinal J (also with reference to the divergence calculation 
inpTcTl). 



1. Longitudinal Equation 



Taking the ux of the \7 x H equation removed a part of the field dynamics lying in the direction of the propagation 
vector u. This part also include the response to longit udinal c urrents. To rectify this omission, we take the dot 
product with u. Using the standard vector identity eqn.( 12.0.13D : 



V ■ (^Ax = A - (y X - B ■ (y X 



and using V x u = 0; we get 

u- (\7 X = --luo? (^u - Ej + u- J 

vector identity: V • x ifj = —iLoa^ * {u ■ Sj + u ■ J 

premultiply: aj.V • (^u x PrH^ = -~iLo(3ra^ * {u- a^E^ + arPrU- J 



arV ■[G+ -G- 



G+ + G- 



2arf3rU ■ J 



(4.2.21) 

(4.2.22) 
(4.2.23) 
(4.2.24) 
(4.2.25) 



This is the same as the source-free divergence- difference equation seen in III C 1 ; where we see that both parts of 
the LHS are zero if u ■ (^G+ + G"^ = and u - J = 0. 



2. Simple Cases 



In the non-magnetic (/3c = 0) case, 



V X G* 
±VG° 



^ILOarPr U X G + lU!ar(3rG° =F — - — ttc * ( W X 
+lLUarf3rUG° 



G+ + G- 



In the time domain, eqn. (4.2.26) is 



V X G± = ±a 



Q-c * /3r 



± MX PrJ 



(4.2.26) 
(4.2.27) 



u X 



G+ +G- 



± ux (3r*J. (4.2.28) 
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Plane Polarized & Magnetic (G^ ^ ± Hy): 

The y components of the LHS and RHS contain the Gx val ues - as a result of the curl and cross-products respectively. 
Compare the previous stand-alone calculation in eqn.(4.1.8) to 



dzG^ = T«war/3r Gt T ''-^Uc * [G+ + G J ^ 



/3. * [G+ - G-] 



(4.2.29) 



This contains a single sign discrepancy compared to eqn (4.1.8); the RHS has the opposite sign. 
NB: recheck curl/cross-product component signs? (checked OK 20041028) 
Vector Moving Frame, Transverse, Magnetic: 



It is worth transforming to the moving frame at this end point, rather than at the beginning as in subsection IV D 



Since eqn. (3. 8. 5) tells us that VQ = V'Q — a//3/u x d[Q, with ^ = aj^f /a^Pr, we get 



V' X G^ -f lUOafPf U X G^ = ^lUJUrPr u X G^ 



ac* X 



G+ + G- 



/3c * Ut X 



G+ -G- 



V' X G± = T«war/3r (1 T u x G± 

G+ + G- 



=P Qfc * I U X 



G+ -G- 



(4.2.30) 
(4.2.31) 



This agrees with the standalone vector moving frame calculation in the transverse case (see subsection IV D); which 
culminates in eqn (4.4.22). I use the transverse case (not the general) to avoid the moving-frame altered V that occurs 
in one of the terms. 

Plane Polarized Moving Frame, Magnetic: 

This is just a special case of the above vector moving frame equation - 



d..Gt = Tiu^arPr (1 T Gt T * [Gt + G"] - * [G+ ~ G"] 



(4.2.32) 



This has sign conflicts with the stand-alone calculation in subsection IV C| which culminates in eqn ( 4.3.14 ). 
Beltrami-like form: 
This might be written 

{{Vx)±iujarf3r{ux)}G'^ = W^, (4.2.33) 
where is a "source" term, but note that actually it will be some complicated function of G^. 

3. Comments 

Note the uiarPr = ^/cf prefactors on the RHS's - these trivially convert {uj/cr — > uj/ujk^^ —> k^) to a wavevector 
prefactor kr, as would be expected for this kind of spatially propagating description. 

Inspection of the equations shows that even if (e.g.) G+ starts identically equal to zero, that it will pick up some 
value due to the interaction of the G~ components with the non- vacuum part of the permittivity. However, so as long 
as Uc <SC OLr, the transfer will be slow, and because of the leading 2ia; term, G"*" will rapidly oscillate. As a result, we 
can assume (in an appropriate limit) that any G'^ contribution to the evolution of G~ will spatially average to zero 
in a kind of rotating wave approximation (RWA), although we should properly define criteria for "slow transfer" and 
"RWA" 

One approach that would avoid throwing away the backward propagating terms would be to describe the G^ field 
with the reverse frame transformation. Whilst this would result in some inconvenience, when simulating numerically 
it could be useful in that fields propagating in the unimportant direction could be thrown away only once they had 
ceased overlapping with the range covered by the pulse travelling in the important direction. 



C. Plane-polarized moving frame (Example only) 



NOTE: In this example calculation, the convolutions required for the nonlinearity are not explicitly included in the 
notation. This does not affect the result, as long as they are reinserted correctly in the final equations. 
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Here I use the (non relativistic) frame translation defined by eqn. (3.8.4) in subsection [II H. Starting with the 
source-free Maxwells' equations, we have - 



dz 
dE,{t) 



d 



dt 



dz 





' d 












' d 











H,, 



Ex — — 



dt' 
d 



at 

--^[^i^Hy+afPfEx] 



(4.3.1) 
(4.3.2) 
(4.3.3) 
(4.3.4) 
(4.3.5) 
(4.3.6) 



Now, because the convolutions e * E^ and ^ * Hy complicate the analysis, I will fourier transform into the frequency 
domain where the linear response becomes a simple product. Using the definitions in subsection ||, I can proceed 
with the calculation - 



d ~ 



dz' 



-Hy{uj) = luj a{ijjYEx{oj) + af{uj)(3f{ijj) Hy{uj) 
P''Hy{uo) + af{u)l3f{uo) Ex{uo) 



^-Ex{^ 



d 



lUJ 



dz' 



■Ex — 



— (3rHy 



dz 



■UrEx = 



lUjarf3r 
lUjarf3r 



arEx + ^PrHy 



+ lUlPrCtc-drEx, 



PrHy + S,arEx + lUJarPc-PrHy 



(4.3.7) 
(4.3.8) 
(4.3.9) 
(4.3.10) 
(4.3.11) 
(4.3.12) 



Adding and subtracting eqns. (4.3.11, 4.3.12), and using the fourier transform of the definitions for = arEx ± 
PrHy, we get 



-r— -G* = ± (1 ± £) ILUardrG^ + lUjar0c-0rHy ± ILOfSrac-arEx 

az' 



± (1 ± lUJarl3r& 



lUI- 



(g+-G-) ±z^^(g+ + G-) 



The "non-magnetic" limit, where /3i = is 

d 



(g+ + g-) 



(4.3.13) 
(4.3.14) 

(4.3.15) 



The reference limit of the above equations (with ai = 0; /3i = 0) when the frame is matched to the reference Cf = Cr 
is exactly as you would expect, given that the new frame is moving in the direction of the G~ component, and since 
the pulse does not evolve (relative to its centre) in a vacuum, dz'G^ is zero (see eqn. ( 4. 3. 17 )). Conversely, the G"*" 
component is moving the opposite way to our frame a t c, thus relative to the frame it moves at 2c, hence the derivative 
is double that in a stationary frame (see eqn. ( [4.3.16 )). 



Ag+ = ±2^G+ 



dz 
^ G- 



(4.3.16) 
(4.3.17) 
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D. Vector moving frame (Example only) 



NOTE: In this example calculation, the nonlinear convolutions are not explicitly included in the notation. This does 
not affect the result, as long as they are reinserted correctly in the final equations. 

Note that is better to follow the (authoritative) vector stationary frame derviation (see IVB) instead; whe re the 
frame translation is applied to the wave equation. I use the (non relativistic) frame translation defined by eqn.(3.8.5) 
in subsection III H ; and using the definitions in section ^ I can proceed with the calculation - 



into moving frame 
into temporal-frequency space — > 

swap transverse — > 

premultiply — > 



u X 



V X H{t) 

V X E{t) 

V X H 

V X E 
V X H{uj) 

V X E 
(V X H^j 

V X E 



sum & diff 



ux [W' X PrH 

V X UrE 

V' X arE ± ux (w' X PrH 
collect terms — > 



+dte*E{t), 
-dtfi * H{t) 

+dt'e * E — oifPf u X dt'H, 
—dffJ, * H — oifPf u X dt'E 

a{ujfE{Lu) - af{uj)/3f{Lu) u x H{uj) 

(3{ujfH{Lj) + af{uj)Pf{Lj) u X E{uj) 

a'^ u X E — afPf u X uxH 

[0^H + afPf uxE 

' PrO^ U X E — afPfPr U X 
^ arP^H + Q!//3/Q!r U X E 

'^aj.0^H + afPfUr u x E 
=F luj ^/3,.Q!^ u X E ~ afPfPr u x 



il X H 



ux H 



^lUJ {o? =F OifPfUr) UX E 



ux H 



The vector fields (defined in eqn ( 3.3.1 )) are 



CirE + U X (3j.H 



(4.4.1) 
(4.4.2) 

(4.4.3) 

(4.4.4) 

(4.4.5) 

(4.4.6) 
(4.4.7) 

(4.4.8) 



This means I need to convert both the second term on the LHS of the sum-and-difference equation above, as well 
as the RHS. It is most imortant for the LHS to be simple, because this will define the type of propagation specified 
by the RHS. To convert the LHS, we need two vector identities, as used in subsection XH to get eqn. ( |12.0.iq ). So 
using 



U X 



(y X -V (u-H^ = \/ X (ux 



(4.4.9) 



WARNING! These are moving frame V's arriving from the vector identities!! 

If I retain the u ■ /3rE term, which would be zero for strictly transverse E fields (TEF), then 



V X arE ± 



u X 



a /3r T CtfPfC^r) U X E 

VxarE ± V X lux f3rH] ±V' (u- f3rH 



U X 



(C) 



V X 



arE ± 



(u X Pr 



H 



^lUJ (oP'Pr T OlfPfar) U X E 

+ ILO (ar P^H ±afPfPr U X 
^lUJ (oP'Pr T CtfPfCfr) U X E 

+ tuj (arP^H ± afPfPr u X 
T V • PrH^ 



ux H 



ux H 



ux H 



(^.4.10) 



(4.4.11) 



(4.4.12) 
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Now note that eqn.( 12.0.12| ) means that 



it X H 



u- H 



u-H 



and so 



X G — u X UrE ±ux 



U X [3rH 



= U X UrE =F (3rH ± U ■ PrH 



(4.4.13) 
(4.4.14) 
(4.4.15) 



So, 



f3^H ±aff3ff3r 



(4.4.16) 



u-H 



^ILO {a^ Pr T OifPfar) U X E — VjJ T CtfPfPr) u X 



u — arfi^u X 



ux H 



u X H 



± afPfPr u X 



ux H 



U-H 



(4.4.17) 

u T V (M-/3ri?)(4.4.18) 



I now separate the interaction parts (depending on etc, /3c) from the reference parts (depending on Or, /3r), an d the n 
substitute (as far as possible) expressions containing G rather than E or H, by referring to eqns ( ^.3.3| ) and (|3.3.4|) . 
Hence 



V X = ^ILOarPr (1 T ^ ^ ^rE — lUOarPr (1 T ^ 

=F tUacPrU X UrE — ILOarPcU X U X PrH 
— ^ILOarPr ['^T C) (^rE ±U X UX^rH j 

=F ILOacPrU X UrE — ILOarPcU X U X P^H 



U X 

tLUarPc 



ILUarPc 



U- Prli 

U- PrH 
U- PrH 



— =Fzwar/3r (1 T C) ^ +lL0arl3r U-(3rH 



u tV 



=F lUjacPrU X arE — tuarPcU X ux PrEl + lUjarPc u - PrH u 



U- (3rH U T V [u- Prii^ 

(4.4.19) 

U ^V'{u-PrH) 

(4.4.20) 
(4.4.21) 



Note that the non-transverse corrections depend only on the Magnetic field H . I assume this is a curiousity of the 
chosen G definition, as (I think) choosing G^ ^ u x E ± H would produce non-transverse corrections depending on 
the electric field. As yet I do not have a simple vectorised form for calculating H purely in terms of G's without these 
non-transverse corrections occuring. 

So using E = {G+ + G-)/2ar] ux H = {G+ - G-)/2Pr] H = u x {u x H) - {u - H)u 



V' X G^ = ^lUJUrPr (1 T " X G^" 



tOJarPr 



U- PrH 



u tV 



'(u-prH) 



lUjacPr 

T — 7: — u X 



G+ + G- 



ILOUrPc 



UX [G+ -G- 



non-magnetic - 
transverse - 



^^ILUttrPr (1 T C) W X =F 
^ILJUrPr (1 T C) W X =F 



ILOacPr 



lUiacPr 



G+ + G- 
G+ + G- 



ILjarPr 



U-PrH 



U - PrH 



(4.4.22) 



u ^V[u-PrH) (4.4.23) 
(4.4.24) 
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E. Simulations: a step by step guide 

1. The medium 

Consider a plane polarized field propagating through a non-magnetic medium. I choose a stationary frame, and 
a dispersionless linear reference (so Ur^Pr are constants). The permeability of the medium is the vacuum value /io, 
and the permittivity has three contributions: a vacuum (reference) eo, a linear dispersion enit), and an instantaneous 
nonlinearity l ■ The medium properties are therefore broken down in the following fashion - 

e = eo + ?f(c^) + ef^ (4.5.1) 

e = al + l^ar + e^^OLr (4.5.2) 

so ef = ef(t^)/ar (4.5.3) 

ef^ = (4.5.4) 

We can see here that choosing a dispersive reference medium would modify both the dispersive correction and the 
nonlinear term. Im the nonlinear case, this adds a time dependence which did not previously exist, which would 
increase the computational effort needed if the nonlinearity was instantaneous. 

2. The full wave equation 



The first order evolution equation, using the parameter definitions at eqns.( 2.1.(: , 2.1.8 )), is 



d.Gt = Ti^arPr (1 T£,)Gt T [Ct + G'] (4.5.5) 

= T^u^arPr (1 T T '-^^ [Gt + G'] T ^ * [Gt + G'] (4.5.6) 

This has three terms - a reference carrier-oscillation- like term ctr), a linear dispersion term ct^), and a 
nonlinear polarization term (~ 5^^)- The straightforward way of solving this using a split step method, where each 
contribution is considered to be relatively weak, hence we solve the above in three steps. Each one integrates forward 
a distance Sz, but since we do one after the other (rather than doing them simultaneously), the procedure is only 
accurate to first order ~ but for small enough 5z, the error can be assured negligible. 

3. The reference evolution 

The reference term just applies a complex rotation to the field in frequency space, which can be calculated exactly 

G^^ = G^(z)xexp TiujarPrilTO ■ (4.5.7) 

This will have no effect if we choose the frame velocity to be the same as the phase velocity of the reference medium. 
This might seem to differ from the usual E field approaches, where choosing the frame velocity to match the group 
velocity gives the best cancellation of propagation terms. However, remember that the way we have chosen our 
reference medium (in this example) means that group velocity corrections appear in the second RHS term as part of 
Off . Choosing a frame velocity equal to the group velocity here will leave a residual reference- like term, but this can 
then cancel with part of the group- velocity-like contribution from the dispersion term, which in most cases would give 



the best cancellation - just as in the usual E field approaches. If we wanted, we could easily rearrange eqn.(4.5.7) to 
incorporate such a cancellation, and then solve the equation appropriately. 
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4- The dispersive correction 



The next step is solving for the linear dispersion — /ar in the frequency domain. Fortunately, this part of the 
equation is (also) trivial to solve exactly in the frequency domain, and results in the following exponential solution - 



G^{z) X exp 



=Fzfc G^.Sz 




(4.5.. 



Although both reference and dispersion steps can be solved exactly using exponentials, there is an important 
difference. The reference evolution of G"*" depends only on G"*", whereas the dispersion evolution depends on the sum 
G+ + G^, since the dispersion acts on the electric field. In the forward-only approximation, of course, G~ = and 
the two steps can be trivally combined - as is automagically done in most approaches solving for the propagation of 
optical pulses. 



5. The nonlinear correction 



The final (third) step is to transform into the time domain and solve for the n-th order nonlinear effects. Since our 
reference ar,(3r are constants, a^^ — e^'^/fir is trivial to calculate and will be just a^^ — x*-"^£^"~^/Q!r, so (for a 
simple Euler method integration) 

Gt{z + 6z) ^ G^j,{z + 6z) ±£^^^_[G+^(z + fc) + G-^(z + fc)]"<5z (4.5.9) 

For a narrow-band field, the time derivative would be dominated by (proportional to) its centre frequency; and 
indeed in most envelope theories we only see a factor of (carrier frequency) in the analogous expression. For a 
wider-band field, there would be corrections, as found for the GFEA approach |Q. 

In an envelope theory, the time derivative would be dominated by the carrier frequency coq, hence we would only 
see a factor of ujq rather than the differential here. A few-cycle envelope theory would give us corrections to the loq 
prefactor. If you wanted to, the derivative used in this nonlinear step could be calculated in the spectral domain - 
you'd have to FT G into the time domain, raise it to the n-th power, and FT back... 



6. Initial conditions: matching a pulse to the medium 



Here I describe how to make the "best matched" initial conditions for G^ describing a pulse propagating only 



in the direction. Note that it may be helpful to read section IX if you are confused by the use of terms like 
"pseudo-reflection" or by the meaning of the co- propagating G~ component. 
The procedure is: 

(1) Pick a suitable electric field E{uj) 

(2) Find out the full dispersive properties of the medium, as described by and fii. You can even put the nonlinear 
or any other properities into and fii as well - but only if you will still know how get a solution for steps (3) and (5) 
below with those added complications. 

(3) Calculate the H{u!) corresponding to E{u!) for a "-(-" propagating pulse: 



= y/7~p)E{Lu) + ^fH{u)H{Lu) (4.5.10) 
Hie.) = J^Eiu;). (4.5.11) 



(4) Choose our reference medium parameters er{uj) and Urito). 
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(5) Calculate our initial for the chosen reference medium, given our initial E{iLj) and H{uj) fields. 

= y^er{uj)E{uj) ± ^/ ^ir{uj)H{Lj) (4.5.12) 



V Mi(^) 



/ 7 — \ / 7 — \" / 

V M^) 



(4.5.13) 

E{uj). (4.5.14) 



These G^ initial conditions are the best guess we can make for a "+" propagating pulse in our medium. 

Note that in step (3) above, I have (in effect) used the definition of G^ with a "reference" equal to the full dispersive 
properties of the medium, and achieved a "+" propagating pulse by setting G~ to zero. In the case where we intend 
to use these full dispersive properties as or simulation reference, of course step (5) becomes trivial. 

Because our reference er,fJ-r will not (in general) include all the dispersion ei(uj), fii{uj), our initial conditions will 
contain both a G^ and a G~ part. This initial G~ part will co-propagate along with the G^, the combination being 
a (hopefully) good match to the medium. 

If our "full" dispersion ti{uj) , ^i{Lo) are not a perfect match to the simulated medium, then there will be an initial 
pseudo-reflection. This pseudo-reflection will consist of a reverse propagating G~ , whose size will be equivalent to that 
of the reflection from an interface between the dispersion used to construct the initial conditions and the simulated 
medium (with + Cc, Mr + Mc)- Mostly this kind of pseudo- reflection can be eliminated with appropriate choice of 
ti{uj), iJ,i{uj); although nonlinear corrections will usually be impossible to treat, and so will still cause pseudo-reflections 
- albiet hopefully very small ones. 

For example, if we construct our initial conditions using a vacum ei{uj) = eQ,fj,i(io) = fiQ, but simulate a medium 
with some refractive index n, this will be equivalent to simulating a ei,fii er,fJ-r interface. Such an interface will 
cause a (pseudo) reflection, which will consist of a reverse propagating G~ , whose size will be equivalent to that of 
the reflection from a vacuum/n interface. A mirror image of the pseudo-reflection will also appear, and co-propagate 
with the pulse - to ensure that the tramsmitted E field has reduced appropriately for such an interface. 

Note that the "forward (co) propagating" G^ component still has a Poynting vector directed backwards. 



7. Comparing alternate references 



Here I do the simplest possible comparison between two different implementations of the same physical system; and 
see that (within the limits set by the approximaions) that they are identical. In the uncoupled stationary frame (G^ 
only) case, the wave equation with only a reference term is 



d^G"^ = =F«wa'o/3r G+ 

= ^iLo [al + araL]^^'^ Pr G^ 

— =F^cJQr/3r [1 + tti/ar]"^^^ G+ 
1 



-aL/ctr 



Gjr 



arf3r 



G ^ 



TiLoarPr G+ T — 7i — 



(4.5 


15) 


(4.5 


16) 


(4.5 


17) 


(4.5 


18) 


(4.5 


19) 


(4.5 


20) 



This is identical to the equivalent wave equation with both a reference and correction term, with the two sets of 
a parameters linked by the appropriate relation e = u'q — + UrUL- The uncoupled approximation (no G~ i.e. 
|G^| <C |G+|) and small correction (1 ^ ai/a,.) are true in the same regime. 
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V. FIRST ORDER EVOLUTION EQUATIONS (PRIMARY H) 



If we wer e more interested in H than E, we might instead define the field variables as in HID, using cqn (3.4.1 
3.4.4 ^.4.5|) . So 



G"±(w) = Priio) H{lu)±ux ar{Lo) E{uj), 

1 r. 



so 



u X arE 



and 



2 

PrH = - 



G'+ + G'- 



(5.0.1) 
(5.0.2) 

(5.0.3) 



A. Vector stationary frame (Authoritative) 

Here I avoid calculations using the x, y, and z components of the field, and retain a fully vector description. Using 
the definitions in subsection I can proceed with the calculation - 



into temporal-frequency space — > 
swap transverse components 

premultiply — 



V X H{t) 

V X E{t) 

V X H{uj) 

V X E{lo) 
V X H 

ux (v X E 



V X f3rH 

u X ( V X arE 



sum- and-difFerence 



V X iSrH ± u X V X arE 



+dte * E{t) + J, 
-dtfi * Hit) 
~VjJ a 

+IUJ (3{ujf * H{uj) 
—luj a'^ * E + J, 
+tLU P"^ * (ux 

-lUJ PrU^ * E + PrJ, 
+IUJ UrP^ * {u X 

— ILU PrU^ * E ± ILO arP^ * 



UX H 



The vector G'^ fields (defined in eqn ( 3.4.1 )) are 

G'± = Prii ±ux arE 



(5.1.1) 
(5.1.2) 
(5.1.3) 

(5.1.4) 
(5.1.6) 



This means I need to convert both the second term on the LHS of the sum- and-difFerence equation above, as well 
as the RHS. It is most imortant for the LHS to be simple, because this will define the type of propagat ion spe cified 
by the RHS. To convert the LHS, we need two vector identities, as used in subsection XH to get eqn.(12.0.1C). So 
using 



ux {\I X H] -\I [u- H 



V X [ux H 



If I retain the u- arD term, which will be zero for strictly transverse D fields, then - 

V X Prli ± ux {v X arEj = -luPra^ * E ± lujarP'^ * {ux itj + PrJ, 
V X PrH ± V X {ilx arEj ±\7 (u- UrEj = -luPra^ *E ± ILOarP^ ^ (^itx + PrJ 

PrH ± (uxarE) = —luPra^^E ±iujarP'^ * [u X H] [u ■ arE 



V X , 

V X G'* = -luj 
Now note that eqn.( 12.0.12| ) means that 



^Prtt^ * i? =F arP^ * {u X itj^ {u ■ arE 



u X 



ux H 



u-H 



u-H 



and so 



uxG' 



U X PrH ±U X 



u X arE 



— UXPrH^UrE ± U ■ UrE U 



(5.1.7) 

(5.1.8) 
(5.1.9) 

+ Pr$.l.l{)) 

+ Pr$.i.n) 

(5.1.12) 

(5.1.13) 
(5.1.14) 
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So, 



V X G = -luj ■! Pra~ * 

2 



-ILO 



{/3,c 



u-E 



* u X 



u X E 



u X H 



} 



U-E 



u =p V ( u • arE 



/3^K1.16) 



I now separate the interaction parts (depending on ac, f3c) from the reference parts (depending on g^ , 0r) i and 
t hen s ubstitute (as far as possible) expressions containing G'^ rather than E or H, by referring to eqns ( 3.4.4 ) and 
(|a|). Hence 



V X G'± 



M X arE 



+lU!ar(3r |m 

+ lOJ^rCtc * [U X 
±lU}ar(3rU X G'^ — 

H -z — etc * w X 



± 



U X 



u ■ arE 



u =fV 



• arE^ 



u X arE 



± ILjarPc * U X f3rH — lLjj3rCi. 



u ■ arE 



U +(3rJ (5.1.17) 



ILOarPrU G' =F VG' 
iLja> 



G'+ + G'- 



± 



r I 



G'+ - G'- 



lUjPrCtc* {u G'°) +^^1.18) 



Here the transverse (TEF) and longitudinal parts decouple, since G'^ is guaranteed electrically transverse, and 
uG'° is guaranteed electrically longitudinal. In any case, the two decoupled equations are - 



V X G'=^ = ±lLUarPr u X G'* + 



-ar * u X 



G'+ + G' 



lUJar . 
± Pc * I U X 



G'+ - G'- 



±VG'° = —ILOarPrU G'° — ILOarPc * (u G'°) . 



+ (irJ (5.1.19) 
(5.1.20) 



We see from the following V A 1 that these two equations are sufficient to describe the fields as long as there is no 



longitudinal H component (i.e. G'^ is transverse); (also with reference to the divergence calculation in HID 1) 



1. Longitudinal Equation 



Taking the ux of the \/ x E equation removed a part of the field dynamics lying in the direction of the propagation 
vector u. To rectify this omission, we take the dot product with u. Using the standard vector identity eqn.( 12^0.131 ): 



V • 



(Ax b) 



= A- V X B] - B ■ V X A 



and using V x u = 0; we get 



vector identity: 
premultiply: 



u- [V X Ej = -zujP^ * (^m • iJ j 
V ■ (ux eJ ^ -iuj0^ *(u-H^ 

PrV ■ {u X arEj = —lUOarP^ * {u ■ PrH^ 



arW ■ (G'^ 



-G' 



= —luia. 



G'+ + G'- 



(5.1.21) 

(5.1.22) 
(5.1.23) 
(5.1.24) 
(5.1.25) 



This is the same as the divergence-difference equation seen in pi D 1 ; where we see that both parts of the LHS are 
zero if w- (g'+ + G'- 



0. 



2. Simple Cases 



In the non-magnetic (/3c = 0) case. 



V X G'=^ = ±iujarl3r uxG'^ + 



±VG'' 



-ILOarPrU G' 



lUjPr 



-ac * u X 



G'+ + G'- 



(5.1.26) 
(5.1.27) 
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In the time domain, eqn.(5.1.26) is 



f3r 



* ar * \ u X 



G+ + G- 



± ux(3r*J- (5.1.28) 



Plane Polarized & Magnetic (G^ ^ ± Hy): 

The y components of the LHS and RHS contain the Gx val ues - as a result of the curl and cross-products respectively. 
Compare the previous stand-alone calculation in eqn.(4.1.8) to 



lLjj3r 

~2 



(5.1.29) 



Vector Moving Frame, Transverse, Magnetic: 

It is worth transforming to the moving frame at this end point, rather than at the beginning. Since eqn.(3.8.5) tells 
us that VQ = V'Q — oirPrU x ^'tQ^ with ^ = arj3r/ UrPn we get 



V' X G'^ + ILOarPr U X G'^ — ±lLUarl3r u X -I- ^ ac * (u X 



G'+ + G'- 



V xG^ ^ ±iuar(3r (1 T uxG 



-etc * W X 



G'+ + G'- 



U X 



G'+ - G'- 



/3c * u X 



G 



1+ 



Plane Polarized Moving Frame, Magnetic: 

This is just a special case of the above vector moving frame equation 



G^ -^r G J, 



G^ 



(5.1.32) 
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VI. FIRST ORDER EVOLUTION EQUATIONS (PRIMARY D) 



Fol lows the de r ivatio n in the previous section, but uses the "Primary D" (or F^) form, as defined in IIIC, with 
eqn. ( 3.5.1 , 3.5.S , 3.5.5 ). So (dropping the tildes out of laziness) 



1 



so 



and 



MX f]-^ B = 



+F- 
F+ -P- 



(6.0.1) 
(6.0.2) 

(6.0.3) 



We would do this so that the displacement field (and magnetic fields) could be easily reconstructed from ; and, 
in analogy to the THF and TEF mentioned above, we would use a transverse magnetic induction B approximation 
(TBF). This this form differs from the primary E form that defines in that it contains a~^D = a~^a'^E not 
aj,E; i.e. the electric-like contribution differs only by e/e^ — c? joi^. A similar comparison holds for the magnetic 
components. 



A. Vector stationary frame (Authoritative) 

Here I avoid calculations using the x, ?/, and z components of the field, and retain a fully vector description. Using 
the definitions in subsection ^ I can proceed with the calculation - 



into temporal-frequency space 



swap transverse components 
premultiply — > 

simplify — 



V X U(t) 

V X E{t) 

V X iJ(cj) 

V X E(uj) 
ux (v X H 



V X E 

M X ( V X Q!r/3^ * H 



+dtD{t) + J, 
-dtB{t) 

—VjJ D{uj) + J(uj), 

+ILO B{lo) 

— ILOU X D + U X J , 

+IUJ B 



^ = —luj arP^ * {ux + UX arP^ 



V X f3ra * E = +IUJ (3ra * B 



(6.1.1) 
(6.1.2) 

(6.1.3) 

(6.1.4) 



sum-and-difference 



u X X arB 

V X 

V X /3rD ± ux (v X arBj 



—luj arl3^ * \ ux Dj + ux ar0^ * J, 
+iuo(ira^^B (6.1.5) 

+IU} PrCi^ * B ^ lUJ ttrP^ * {u X 

± uxar0^*J, (6.1.6) 



The vector F^ fields (defined in eqn (3.5.1)) are 

F^ = a-^D±ux (3-^B 



(6.1.7) 



This means I need to convert both the second term on the LHS of the sum-and-difference equation above, as well 
as the RHS. It is most imortant for the LHS to be simple, because this will define the type of propagat ion spe cified 
by the RHS. To convert the LHS, we need two vector identities, as used in subsection XH to get eqn.( 12.0.1C ). So 
using 



u X 



(v X -V (^u- = V X (^M X i?) 



(6.1.8) 
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If I retain the u x [3rH term in the above, which will be zero for strictly transverse fields, then - 

/3^*J (6.1.9) 
ux 13^* J (6.1.10) 
(ux =F V ( 

±uxarP'^*J (6.1.11) 
luj^PrU^^B T arP'^*{uxD^^ T^{u-arB^ ± U X UrP^ * (B.1.12) 



V X PrD ± UX (V X UrB 
V X PrD ± V X {uX UrBj ±V {u ■ ttrBj 

PrD ± {ux arB 



V X 



V X F± 

Now note that eqn.( 12.0.12|) means that 



ux H 



it- H 



it- H 



and so 



(arPr) U X = M X PrD ± it X 



U X UrB 



U X PrD q= UrB ± 



U • UrB 



(6.1.13) 
(6.1.14) 
(6.1.15) 



So, 



{arPr) VXF^ 



ILU < Pra * 



v.- B 



U — PrO? * it X 



it X B 



± it X arP'^ * J 
^UrP^ H X D ± Pra^U X 

± U X Q!,./3^ * J 



Ux B 



arP"^ * (itx D^^ 

ILoPrO^ 



TV 



{u ■ UrB^ 



it- B 



it =fV 



• OrB^ 



(6.1.16) 



(6.1.17) 



I now separate the interaction parts (depending on etc, Pc) from the reference parts (depending o n g^, Pr), an d the n 
substitute (as far as possible) expressions containing F rather than D or B, by referring to eqns ( 3.5.3| ) and (3.5.5). 
Hence 



V X 



^rPr I 



it X a^^D ±it X 



ux p-^B 



IbjarPr 



u- p-^B 



Pc* {ux ^D^ — iLuPrac * x u x P^ ^B 

± UXp-^{Pl+PrPc)*J 

^ILOUrPr U X F"^ + ILUarPrU F° =F VF° 

ILUP, 



it tV 

imPrOLc * 



B 



u-p;^B 



(6.1.18) 



T -T^Pc * M X 



F+ + F- 



-Ur * \ U X 



^ + lUlPrOic * {U F°) 



± U.x{Pr+ Pc)*J 



(6.1.19) 



Here the transverse (TDF) and longitudinal parts decouple, since F^ is guaranteed displacment transverse, and 
ilF° is displacment longitudinal. Note that in this transverse/ longitudinal split, the role of the current J has not 
been checked. In any case, the two decoupled equations are - 



V X F — ^lujttrPr it X F =F — - — Pc * yit X 
±VF° = +iujarPrit F° +iujPrac* {it F°) . 



F+ +F- 



ILoPr 



-OLc * U X 



F+ - F- 



± Ux [Pr + Pc)* J 



We see from the following VI A 1 that these two equations are sufficient to describe the fields as long as there is no 
lo ngitudi nal D component (i.e. F^ is transverse) or longitudinal J (also with reference to the divergence calculation 
in llHEll). 



1. Longitudinal Equation 

Taking the Ux of the V x H equation removed a part of the field dynamics lying in the direction of the propagation 
vector U. This part also include the response to longitudinal currents. To rectify this omission, we take the dot 
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product with u. Using the standard vector identity eqn.( 12?0.13D 



V • 



= A- V X B -B- 



(v X A) 



and using V x u = 0; we get 



vector identity: 
premultiply: a~^V • {u x (3^^ B 

V • [f+ - P- 



V ■ (ux H 



— ILUU ■ D + u ■ J 

— ILUU ■ D + u ■ J 



(6.1.22) 

(6.1.23) 
(6.1.24) 
(6.1.25) 



-uvar iPr + /3c) * (^i 



F+ +F- 



+ {(3r+Pc) * (m- j) (6.1.26) 



This is the same as the divergence-difference equation seen in III E 1 ; where we see that both parts of the LHS are 
zero if u ■ (^F+ + F"^ = and u - J = Q. 



2. Simple Cases 



In the non-magnetic (/3c — 0) case, 

V X F"^ — ^lUJarfir u X F^ 



-Uc * \ U X 



F+ -P- 



± U X {Pr + Pc) * J 



±VF° = +tLuarf3rU F° +iujPrac*{u F°). 



(6.1.27) 
(6.1.28) 



In the time domain, eqn.( 6.1.27 ) is 

V X = ±dt {ar * Pr) * (u X P^^ 



* Q!r * U X 



F+ -F- 



± u X (/3r + /3c) * J.(6.1.29) 



Plane Polarized & Magnetic {F^ r^.j, ±By): 

The y components of the LHS and RHS contain the Fx val ues - as a result of the curl and cross-products respectively. 
Compare the previous stand-alone calculation in eqn.(LL^) to 



dzF^ = ^ILUarPr F^ 



-^a, * [F+ - F-] T ^/3c * [F+ + F'] 



(6.1.30) 



Vector Moving Frame, Transverse, Magnetic: 

It is worth transforming to the moving frame at this end point, rather than at the beginning. Since eqn.(3.8.5) tells 

us that VQ = V'Q — a^PrU x d'lQ, with ^ = ar(3r/ar(3r, we get 



V' X F -\- liOUrPr U X F ~ ^ILUarPr U X P =F — —^f3c*[UX 



V' X F± = =F«wa^/3r (ItC) uxP"^ 



-/3c * U? X 



F- 



p+ + p- 



-Ur * M X 



p+ - p- 



-ar * U X 



F+ -F- 



(6.1.31) 
(6.1.32) 



Plane Polarized Moving Frame, Magnetic: 

This is just a special case of the above vector moving frame equation - 

d.,F^ = T^a.(3r{lT0 ^ 0^ , {[F+ + F']) - '^a^ , ([f+ - F']) 

Beltrami-like form: 

This might be written 

{{Vx)±tujar(3r{ux)}P^ = W^, 
where W"^ is a "source" term, but note that actually it will be some complicated function of G^. 



(6.1.33) 



(6.1.34) 



VI. A 



32 



1^1 



DETAILED CALCULATION 



Directional fields: 



Dr.Paul.KinslerQphysics.org 

http: / /www.kinslcr.org/physics"^ 



VII. FIRST ORDER EVOLUTION EQUATIONS (PRIMARY B) 



Foll ows the derivation in the previous sections, but uses the "Primary B" (or F'^) form, as defined in IIIF, with 
eqns. (l3.6.l| , ^.6.4 |3.6.5|) . 

NOTE: Not yet calculated, but it will have many similarities to the preceeding three derivations in sections |fv| 0, 

v3. 



VIII. THE ENERGY DENSITY AND THE POYNTING VECTOR 



Here I discuss the role and calculation of the energy densities and Poynting vector for the EM field described in 
terms of the Primary E (G^) fields. While it turns out that the THE Poynting vector is always a simple function of 
the G± fields, the energy density is more complicated and depends on the degree of mismatch between the reference 
parameters and the true material parameters. Because of this difhculty of reconstructing H f rom G^, all these 
calculations currently resort to restricting themselves to the THE case only (i.e. G° = 0, see [II C ). 

In practise it might well be simpler to just reconstuct the E and B fields from known G^, and use the traditional 
forms for energy density and Poynting vector rather than evaluate the complicated expressions derived in this section. 
Nevertheless, these complicated expression may provide insight into the nature of the G^ fields, even if they a not 
used for computation. 

Here I will restrict mself to the non-dispersive reference case, where e, are constants independent of frequency. 
Without this, getting E{t) from G^{t) requires a deconvolution to get the frequency-dependance of e, ^ out of G^{t). 
This restriction can be avoided entirely by taking the formulae below to hold in the frequency (rather than time) 
domain. 

I define 0, denoting a convolving dot product, and ®^ denoting a convolving cross product, so 

(8.0.1) 
(8.0.2) 

I've used 



AqB ^ J duj' A{uj') ■ B{uj ~ uj'); 
A®B = / duj' A{lu') X B{lo -Lu'). 



A{t) * Bit) ■ C{t) ^ A{uj)B{uj) G(w) 
A{t) * B{t)\ X C{t) <-> \a{uj)B{uj) \ C{uo). 



To avoid cluttered equations, I will also use this shorthand notation 



G 



G±/A 



1/2 



(8.0.3) 
(8.0.4) 

(8.0.5) 



A. Energy Density 



The standard expression for bulk energy density of the EM field in terms of the electric and magnetic fields can be 
easily converted into one in terms of G^ in the THE case. In the following, I explicitly allow for the possibility of 
mismatched reference and material parameters (i.e. e 7^ etc). 
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1. Miscellaneous Identities 

Where, at the transverse approximation: 

{uxA^-[uy.B^ =A'-{uy.lPj = - A ■ (^B y. =-B-[uy.A'^ (8.1.1) 

= -B ■ [uy^ uyc = -B ■ (ux ux (8.1.2) 

X A 



transverse. 



= +B- 
= +B-A- 
(^ux (^ux B^ w +B- A 



B ■ u 



X 

u ■ A 
A-B 



= +B-[[u-u]A- u-Auj (8.1.3) 

(8.1.4) 
(8.1.5) 



The following brief calculation illustrates the non-trivial nature of getting a simple expression in the frequency 
domain - 



Ua = €*E-E 

= J e{t')E{t - t')dt' 



E{t) 



e{t')E{t - t')dt 
= J e(<') J e-"^'E{t - t') ■ E{t)dt 
= J e{t')F{LO,t')dt' 
where F{uj,t')^ e-"^*E{t - 1') ■ E{t)dt 



E{t)dt 
dt' 



(8.1.6) 
(8.1.7) 

(8.1.8) 

(8.1.9) 

(8.1.10) 

(8.1.11) 
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2. Transverse Case 



The time domain energy density calculation is (making use of eqn.(B.0.5)) 



freq. domain 



transverse approx eqn( 8.1.5 ) 



time domain. 



U ^ -e* E ■ E H ■ H 

2 2^ 



2 Uiu) 
8 U{uj) 
8 U{uj) 



eE]QE+ (nH) QH 



eG 



Gt + G 

0G+ + 







Gp + Gf 



+ u X 



eG: 



)o(. 



Gfj. + G^ 



.1.12) 
.1.13) 
.1.14) 



G7 + eG7 G+ + eG7 G, 



Uiio) = 



^G 



■ q(ux G+) -(ux ^G+) X g;) 

u X ^G;) (m X G+) + X /iG;) (u X G; 
G+ + (eG+) G7 + (eG;) G+ + (eG^) G7 
mG+) G+ - (AiG+) g; 

- (a^g;) 0G+ + (mG;) 0G; 

G+ + (mG+) G+ + (eG7) G7 + (mG;) G; 
eG+) G: - (mG+) QG- + (eG:) G+ - (a^G;) 





G^ 







±G- 



^G+ 



^G- 







-G- 



ar 



A 



-G" 







-G+ 







-G" 



■G^ 







8 Uit) 



* G^ 



^G- 





^G+ 

ar 






G- 


— *G+ 


+ 






■ 1 



■G" 

' 1 



*G^ 



G- 



— *G" 



— *G- 



— *G+ 



— *G" 



— *G" 



— *G^ 

a,. 



■^*G- 



The first two lines of eqn.( 8.1.19 ) are the G+^ and G~^ terms that might naturally be expected to occur; the 
final two lin es are G"*" • G~ terms caused by interference between the two fields. In the dispersionless reference case, 
eqn.(B.1.19) reduces to 



freq. domain. 



time domain. 



U(tj) = 



8lt(i) = 



( 




G+^ G+ + 

















67- l-if 

Mr 



G+ I 0G 
* G+ 1 • G 

G+ 1 • G" 



It. 

Mr 



G- G^ 



* G" 1 • G" 

* G" 1 • G^ 



Mr 



.1.20) 



.1.21) 



.1.15) 



.1.16) 



.1.17) 



.1.18) 



.1.19) 
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In the fully dispersionless case, eqn.(B.1.19) reduces to 
freq. domain. 



time domain. 
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.1.221 



.1.23) 



We see here that the energy density takes a more complicated form if we choose to use reference parameters that 
do not match the true medium ones. If the reference parameters have a different degree of mismatch between e and 
fi, which will in fact normally be the case, then the cross terms dependent on G+ • G^ do not vanish. 



3. Full Vector Case 



The energy density in the fully vectorised case can be found by retaining the 



B ■ u 



•■A 



corrections removed by 



the transverse approximation. Thus the non-dispersive reference energy density, as specified by eqn.(8.1.21), becomes 
811 = 
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The fully dispersive version of the corrections is 
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u • I ^ * G^ 



/3r 



*G" 



M • I — * G" 

Pr 



w • I * G- 



^ *G- 



u • ( — * G+ 

Pr 



|(.8.1.24) 



.1.25) 



B. Poynting vector 



The Poynting vector defines the energy flux, and is 

P = E X H 



(8.2.1) 



We will see that a G"*" component always corresponds to a Poynting vector (energy flux) in the -\-z direction, and 
the G~ to an flux in the opposite direction (i.e. —z). This is the basis of the assertion that these G* variable are 
direction, and thus correspond to forward and backward directed field components. 

NOTE: A critique of Poynting vectors - see e.g W. Gough, "Poynting in the wrong direction"^ Eur. J. Phys. 3, 83 
(1982). 
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1. Vector Identities 



Because the convolution integrals are independent of the vector properties, we can work in the transverse limit 
to get 



and 



G+ (g) M X G ^ - 



G~ (g) it X G+ = 



G+0G 
G+0G 
G- QG+ 



G+- 



G- 



G- 



G+ 



= + G-QG+ u 



(8.2.2) 
(8.2.3) 
(8.2.4) 
(8.2.5) 



Transverse Case 



The Poynting vector defines the energy flux. To get a simple expression, we need to work in the transverse case and 
specialise to a non-dispersive reference medium. Some expressions remain in the weaker partially transverse "THE" 

case - see III C , where u ■ H = and 2/?,.iJ = u x ( G+ — G^ ) . Starting from the standard expression, and making 



use of the notation in cqn.(B.0.5)), we get 
frequency domain S{uj) — 



THE approx. 



simplify as below (VIII B 1) 



time domain; 
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+ G, 



Gr g: 



u + O 



uxG+ 



-g: 



u X 



G7 -u 



* g: 



G'^ — G G , 



(a-^G+^ (p;^G+^ - (a-^G-^ (/3,^'G-) 



45 = 



a-^ * G^ 



f3-^ * G+ ) - (a-i * G- ) • * G- 



(8.2.6) 
(8.2.7) 
(8.2.8) 
(8.2.9) 

,10) 
(8.2.11) 
(8.2.12) 
(8.2.13) 
(8.2.14) 
(8.2.15) 
(8.2.16) 
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3. Minkowski Version 



Starting from the standard expression, I use the shorthand notation — G^e.er and G^,j^ ~ G^fi.^r 



avoid cluttered equations; and so 



frequency domain : 



THE approx. 



~4 5* 



M 



simplify as below ( |VIII B l| ) 



time domain; 



to 
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-4.Sm = 



±,G^].(J^,G^]-(±,G-]-(J^.G- 



IX. INTERPRETATIONS 



The most comprehensive discussions on the interpretation of G^ fields are contained in Kinsler et.al. 2005 
Other important remarks on the behaviour of the propagating fields at interfaces are covered in my report "Causality 
in spatially propagated optics" (tb arXiv). 
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X. SECOND ORDER EVOLUTION EQUATIONS 



In sections |^ and ^ I calculated various forms of first-order propagation equation for the field variables. 
However, many pulse propagation theories start from a second-order form for E, as in e.g. M, 



(9f 



Wl)E{r,t)-\df f dt' 



e{t-t')E{f,t') 



(10.0.1) 



To address this point of interest, here I follow a similar path and derive a second-order propagation equation. 
As expected, it looks rather si milar to th e standard equivalent s for E , but with som e ex tra derivative (curl) term s. 
First I attempt to make eqns.( 3.1. (: , 3.1.7 ) look similar to eqn.( lO.O.l ) in subsection X A . Then in subsection XB I 
present a more complete derivation of the vector form in the transverse field case. 



A. Scalar Form 



This is intended as a quick, intuitive calculation to get a scalar second-order equation for fields. It is based 
on the old Fleck||l[ definitions, so I leave the linear dispersion (from the permittivity e or permeability fi) inside the 
polarization P rather than separate it out as an extra term (thus it contains a hidden convolution). This approach 
has been superceeded (I now p ut everything in e, ^), but this part is rath er old and has not been updated accordingly. 
To start, I multiply eqn.(3.1.6) by the LHS derivative term seen in eqn.( |3.1.7| ): 

[dt-d,][dt+d,]G+ = [dt-d,][-dtP^a{G+ +G-)] (10.1.1) 

= - [dt- d,] dtP -<j[dt- d,] G+-a [dt - d,] G- (10.1.2) 

= - [dt- d,] dtP -<j[dt- d,] G+-a [-dtP - aG+ - aG~] (10.1.3) 

= -[dt-d^-<T]dtP-(j[dt-d,-a]G+ +<j^G- (10.1.4) 



(A) 

(B: sub ( |3.1.7| ) for last RHS term) 
(C: iterate) 



(D) 



dtP 



Q = 



N 

dt-d, + Y^ i-T) 

n=l 

Q] dtP -a[dt- d, 
, for (T < 1, 



N 



[dt - d 

AT 

— a 



dt-d, 
Q] G+- 



n=l 



G^ 



a 



(7^+iG(10.1.5) 

(10.1.6) 
(10.1.7) 



It would now be straightforward to introduce a variety of useful approximat ions, e.g. ct <C 1 or similar. Al ternatively, 
I could group the a terms with the derivative operators in eqns.(3.1.6 3.1.7). I rearrange eqns.(3.1.6 3.1.7) to give 



[dt + d, + c7]G+ = -dtP-aG-, 
[dt-d, + a]G- - -dtP-aG+, 

and apply the differential operator from the LHS of eqn. ( |10.1.9| ) to eqn. ( |1 0.1. 

[dt -d, + a] [dt +d,+a]G+ = [dt -d, + a] [-dtP - (tG"] , 



{[dt-d, + a] [dt + d. 



= -[dt-d,+ a] dtP - n [-dtP - cjG+] 
= -[dt-d,]dtP + cT^G+ 
a]-a^]G+ - - [dt-d,]dtP. 



Now I simplify the term in braces {...}, but retain the ordering of cr's and i9, because a contains e 
make the medium dispersive, even though it is customary to add dispersion at a later stage. 



(10.1.8) 
(10.1.9) 

(10.1.10) 
(10.1.11) 
(10.1.12) 
(10.1.13) 

I might like to 





[dt^d,^ 




-a']G+ 


(10.1.14) 




{[{dt + n)-d,] [{dt + c7) + 


dz] 


-a^]G+ 


(10.1.15) 


{{dt- 


f nf + {dt + a)d,- d, {dt + a)- 


dl 


-a^]G+ 


(10.1.16) 


{(an 


r cr)^ + dtd:, - dzdt + adz - d^a - 


dl 


-a^]G+ 


(10.1.17) 




[{dt + <jf + {adz - dza) - 


dl 


-a^]G+ 


(10.1.18) 



[dj + adt + dta + {adz - dza) - dl] G^ 



(10.1.19) 
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And for G , the term would be {d^ + adt + dtcr — {adz — dz<j) — c?^} G . Apply this to the full equation, assert 
the usual case of a uniform material (so ad^G^ ~ dzCrG^), and then duplicate the calculation for ... 



{d^z - dt - ^dt - dta} G^ 
{dl-dl~adt~dta)G- 



= [dt-d,]dtP, 

= [dt + dz]dtP, 



(10.1.20) 
(10.1.21) 



Compare this to the "standard" second order wave equation eqn.(lO.O.l): (1) the d^G^ is analogous to d^E] 
(2) [dt + dz) dtP has an extra dzdtP part which might be expected given the directional nature of G*; (3) the 
— [df + adt + dta) G^ is the analogous term to —df J eEdt; remem ber that a is a condu ctivity term not usually 
included in pulse propagation problems. If we add and subtract eqns. ( |10.1.20| ) and ( |l0.1.2lD , setting a — 0, we get 



G- 
G- 



= d^P, 
- dzdtP 



(10.1.22) 
(10.1.23) 



Of course cx G+ + G , and Hy oc G+ — G , so the equivalence between cqn.( 10.1.22 ) and eqn.(lO.G.l) is to 
be e xpected . Presumably there would also be a related equivalence between eqn.( 10.1.23 ) and the Hy version of 
eqn.(|lO.O.lD. 



B. Vector Form 



This is a vectorised derivation of the second order wave equation for G^ fields. Unlike the previous scalar derivation, 
it is reasonably complete and without stringent approximation (barring its transverse field approx) . The polarization 
term P remains as part of the dispersion term, thus simplifying the equations; until the end when I split it off for 
illustrative purpos es. Fr om section [II , the unit vector in the propagation direction is u, and the vectorised version of 
G is given by eqn.( 3.3.1 ) - 



In frequency space. 



(-) 

-ILU [jl * H 



G^ = drE + UX PrH. 



-lU! ((5: + drCtc) E = W X H - J, 



H = -V X E. 



(10.2.1) 

(10.2.2) 
(10.2.3) 



In the following calculation, I omit the convolution symbols that should appear between the dc, (3c parameters 
and the G^ ± G~ fields; this is in the interests of both my laziness and reducing notational clutter. I split (as 
before in subsection O) the permittivity and permeability into instantaneous and time-dependent parts with e{uj) = 
a{uj)^ = a^ + araciujji&nd similarly for /i and (3). This mimics traditional derivations, where dispersion and nonlinear 

polarization are handled separately. Using u x V x Q — V (u ■ Qj = V x u x Q, I take the scaled time derivative of 



the definition of G^, 

(— iw) arPrG^ = —iLoalPrE ± (— iw) arPl U X H 

= V X (3rH — PrJ + lOJacOlrE =F W X 

= V X PrH — (3rJ + lUJUcarE =F V X 



V X UrE 

u X a-rE 



± lUJU X f3c(3rH 



(10.2.4) 
(10.2.5) 



TV 



U ■ UrE 



= V X (3rH - (3rJ T V X 



U X UrE 



TV 



u ■ arE 



= TV X 
= TV X 



u X arE =F (3rH 



(3rJ TV 



U X arE ±UX U X PrH 



± V X 



U ■ UrE 



(u-PrH) 



± lUJU X PcPrH (10.2.6) 
+ iLuacarE ± luju x (3c(3rH {10.2.7) 
+ lojacarE ± luju x PcPrH (10.2.8) 



PrJ TV 



u ■ arE 



lUJacarE ± lUJU x PcPrH- 



(10.2.9) 
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Now I combine -E's and H's to form G^'s, 

= ipV X uxG^ 
V X G± 



± G°V XT? =F V 



U ■ UrE 



+ lUjacOirE ± ILOU X PcPrH — (3rJ 



-Oic 



G+ + G" 



q=u X 



± V 
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4v 



=Fw X 
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u-[G+ +G- 

ILO 



13 r J 



G+ + G- 



G+ -G" 



2 



-On 



G+ + G- 



2 



G+ -G- 



^t • (g+ + G" 

G+ -G- 



u- (g+ -G-) 



(10.2.10) 

(10.2.12) 
^(3A4-13) 



and apply another scaled time derivative —iLuarPr, 

,2 



u af.p^G'^ = zbiwar/Sr-u X V X G* + —arPrac 



(B) 



(C) 



(D) 



(E) 



(F) 



(G) 



lU) 



u- (g+-G-) 



=FU X V X ( —uvarf^rG ) + — Ckr/Sra 



(=Fu X Vx) =pu X 



u • I G+ - G" 
V X G=^l 



G+ + G- 
G+ + G- 
G+ + G 



T -Y'^rf3rl3c 



G+ -G- 



G+ -G- 



2 



G+ -G" 



M X V X M X 



=F -Vu X V X 
2 



a; 



M X V X U X 

1, 



G+ + G" 

V X G± 

u • (^G+ - G 
G+ + G 

V X G± 

G+ -G 



T yacW X V X 



—arPrac 



G+ -G 
G+ + G 

± U X V X PrJ 

G+ -G- 
G+ + G 



- y/3c?^ X V X 



u • ( G+ - G 

G+ -G 



It • I G+ - G" 
+ iujar(3lJ 



(10.2.14) 
(10.2.15) 

(10.2.16) 



- Y^rPr V 



T y"cu X V X 



u • (^G+ - G- 

G+ + G- 



luoarPlJ (10.2.17) 



ILO 

- -l^PcU X V X 



G+ -G- 



=F -Vu X V X 
2 



M • G+ - G" 



= u X \7 X u X 



± (-iw) ■ 



w X V X u X 



V X G± 



— (— itj) — 
^ ' 2 



— (—lUj) =F -u X Vx 



u • I G+ - G 

G+ + G 



iLoarPrJ ±u X V X /3r^0.2.18) 



(— lu;) =F w X Vx 



{ — ILj) ^ u X Vx 



/3c 



G+-G 
G+-G- 



Cr 



(— zw) =F w X Vx 



f3rJ 



(10.2.19) 



V X G± 



— {—tu;) =F M X Vx 



(— iw) q= u X Vx 

















G+ 


+ G- 


±/3c 


G+ 


- G" 



G+ -G- 



(10.2.20) 
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Time domain (non-dispersive reference) 



= u X V X u X 
at'' 

1 d 
~ 2di 



V X G± 
1 d 



dt 



=F M X Vx 



1 ^ 

=F u X V X 

Cr dt 



G+ + G- 
G+ -G- 



G+ -G- 



f3rj} 



(10.2.21) 



Since I've already assumed transverse fields, in the source-free case V-£' = 0)soVxVx(3^ V^Q. Thus, 



c2 dt^ 2dt\crdt ^ ^ 













1 |ac * 


G+ 


+ G- 


±l3c* 


G+ - G- 



1=0. (10.2.22) 



Note the strong similarity to eqns.(10.1.2C,10.1.21). The usual polarization term P is contained in the (and Pc) 

0) gives us: 

=Fw X Vx Ip. (10.2.23) 



terms. If I write Uc — + oi^^ and P = ara^^E; dropping the magnetic terms (i.e. /3c — 0) gives us 

I d { I d 



c2 dt^ 2dt\crdt ' " 



G+ + G- 



2ar dt 1 Cr dt 



The l/ur prefactor occurs on the RHS because we are relating G^ to the polarization P, rather than the usual E 
and P. 

We could use these second order equations to apply diffraction to a standard first-order equation propagation, by 
split-stepping the diffraction given here V^G* with the usual dispersive and nonlinear propagation. 

C. GFEA Form 

My "Few cycle pulse propagation" detailed calculation ^ has included terms for G^ field envelopes; where the 
GFEA equation is (g = ±1): 

l±uj'aT \±i<7'Ot \ Pr J 2/?^ (1 ± zcr'ar) 



+ — 1 ± @ iga dr®g{l- iqd^)] — — — p-^B{rjL,t t; A) + — — — -— . 

(1 ± ict'Ot) 1 ± la'Or 



(10.3.1) 



Here the @ and @ symbols act like the usual ±,=F ones; but varying according to which G variable is being 
described - i.e. G ® . The normal ±, =F ones are relate to the carrier direction. 



XI. ENVELOPE PROPAGATION EQUATIONS 

It is common in relatively narrowband cases to solve for the propagation of field envelopes centred around chosen 
carrier frequencies rather than the full detailed oscillations of the EM field. In fact, if sufficient care is taken with 
approximations, and if the system simulated is well behaved, even quite wideband pulses can be sucessfuUy modelled 
using envelopes. The advantage of the G^ wave equations presented in sections pv| ^ and VI is that we dispense 
with the second order form of the wave equation, and generate envelope equations from first order equations. This 
greatly reduces the number of approximations without increasing the complexity of the theory. Although a full model 
requires four envelopes to describe the G^, just as a full Maxwell theory requires four envelopes, (both backward and 
forward travelling E and H), this is rarely a case we are interested in solving, and so in practice we can use just one 
envelope. 

A full expansion of G* into forward and backward envelopes G^ , G^ would be 

G±(a>) = Sf{ioTu;o)e^""''- + 9f{uJo±Uj)e^""'' 

+S±(a> T u;o)e^'''"' + g±(^o T u;)e±*'=«^ (11.0.1) 
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where we have suppressed the z argument on the envelope functions for brevity. Note that the argument shift 
(a; ^ a; ± wq) is how the —iloqI part of the carrier wave is accounted for. If inserted into the wave equations, this 
expansion would result in a large number of terms to consider, even for the relatively simple case of a third order 
nonlinearity. However, by specializing to the typical case where we are only interested in forward travelling waves, we 
can instead use 

G±M = g±(c^T^o)e^"=«^ + S^(c^o±^)e^''"^ (11.0.2) 
= g±(c^T^o)e±''=«" + S^(c^T^o)*eT''=«^ (11.0.3) 
G±(t) = g±(i)e±*'=°^e=F*"'«* + g^(t)*e=F*'=''^e±'"°* (11.0.4) 

This has the nice advantage that it eliminates and backward travelling waves, and so we can propagate pulses 
efficiently in a moving frame. This is important, because the backward parts in a moving frame move at twice the 
frame speed. In a full (non-envelope) simulation, we need to either take special efforts to filter out these backward 
parts, or any hoped-for numerical gains are reversed by the finer z-step required to integrate (the backward parts) 
accurately. 



A. Foward travelling envelopes in a first order wave equation 

I start with the plane polarized first order wave eqn. ( 4.2.29|) 



d,,Gt = T^^oc,^5r{\Ti) G± ^ * [G+ + G"] 



/?, * [G+ - G-] 



(11.1.1) 



I insert envelopes using the shorthand notation S — ik^z. Since etc niay contain a dependence on the fields G^, I 
split it using ac = Oc + a*. In the case of dispersion, Cc — ad'^,] but for a nonlinearity it will be the appropriately 
carrier-matched, positive frequency part of Oc, once the field values it contains have been expanded in terms of the 
envelope and carrier functions. I will proceed in the stationary frame case (^ = 0), so 



9. 



e±-a,g±±zfcoe±-g± 



■eT-9,gt *T»fcoe^"gf * = 



lUJ {ttc + a*) Pr 



g^- 



gt *eT= 



fe^^ + g^ *e^- . 



* 



3fn}i.2) 



5j *em 



split c.c. halves; e^^^^/g^ ± zfcoe^^g^ 



TT^ar(3r9fT ^ |g/+g/ I 



cancel e 



dzSf = Tiij^arPr - ko) g^ T 



{g^ + gj*}, 



}l.3) 
(11.1.4) 
(11.1.5) 
(11.1.6) 



If we do the sensible thing and relate the carrier parameters /cq and ujq using the reference parameters (so fco = 
ujo/cr{uJo)), we have wa^/Jr — ko = {to — wo)Q!r/3r + ujoarf3r — ko = {to — uJo)arl3r: 



dzSf = T« (w - Wo) g^ T 



(11.1.7) 



Now transform into the moving frame. Note that we will no longer want the frame to move at the phase velocity 
at ujQ, because that job has been taken over by the carrier; but we could use it to absorb any residual (reference) 
dispersion affecting one of the envelopes (most likely for g^). However, it will then make the residual (reference) 
dispersion term acting on the other envelope (e.g. 5J) twice as large. Anyway 



dz,9f-iu;af(3f9f = Ti{ujarPr-ko) g± ^ ^^^^y^ {g^ + gj *} 



(11.1. 



9,,g± = Ti{ujarPrT^afPf-ko) g± ^ |g± + gj *} (11.1.9) 
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B. Envelopes and the Second Order wave equation 

See which includes correction terms appropriate to fields. Note that there I can think of no reason to use 
the second order form, since the first order form performs as well with fewer complications. 

XII. VECTOR IDENTITIES 

(1) Identity (old Physics 31.320 vector identity sheet, (1-5)): 

v(u-h) ^ {u-V)H + (h -v) u + ux (v X h) + H X {V X u) (12.0.1) 



rearrange u x (v x H^j = V (u ■ H^j - {u ■ \7) H - (^H ■ u - H x {\7 x u) (12.0.2) 

constant unit vector = V (u ■ H^j - {u ■ V) H - - (12.0.3) 

transverse field = - {u ■ V) H - Q - (12.0.4) 

Useful: {u-V)H ^ V (u-h) -ux (v X h) (12.0.5) 



(2) Identity (old Physics 31.320 vector identity sheet, (I-IO)): 

yx(^uxH^ =u (^V-H^-H (V-u)+(i?-v) u~{u-V) H (12.0.6) 

constant unit vector ^ u (^V • J?^ - + - (u • V) H (12.0.7) 

no monopoles = 0-0 + 0- (m-V) H (12.0.8) 

Useful: {u-\/)H = -\/ X (ux h) . (12.0.9) 



These two identities, specialised to the case involving our unit vector u in the propagation direction, mean that we 
can equate 



u X 



(vxH^-v(u-H^ =Vx(uxH^ (12.0.10) 
Note also 

axbxc = {a-c)b- (a-b)c (12.0.11) 



so in our case where a and b are both unit vectors 

ux(^uxX^ = (^u- X]u- {u-u)X (12.0.12) 

Also, I use (old Physics 31.320 vector identity sheet, (1-7)): 

V-(Axb) = A - (y X B] - B ■ {y X A] (12.0.13) 
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XIII. UNITS 



[^] = N/A^ = m.kg.s"^A~2 
[B] = Y.s.m-^ 

[V.A.s2.m-2]2 .m-i.kg^i 

[e] = F/m = A^.s-^.m-^.kg-i 
[E] = Y.m-^ 
[eE^] = A2.s4.ni-3.kg-i .V^.m-^ 

[V.A.s^.m-^]^ .m^^kg"^ 



(13.0.1) 
(13.0.2) 

(13.0.3) 
(13.0.4) 

(13.0.5) 
(13.0.6) 
(13.0.7) 

(13.0.8) 
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